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MATHEMATICAL MODEL OF POPULATION INTERACTIONS 
WITH DISPERSAL: STABILITY OF TWO HABITATS WITH A PREDATOR 


H. Et-Owaipy AND A. A. AMMAR 


Department of Mathematics, Faculty of Science, Al-Azhar University, Nasr City 
Cairo, Egypt 


(Received 21 April 1987; after revision 26 July 1988) 


A system of differential equations of dispersion between two populations 
in habitats separated by a barrier with a predator feeding indiscriminately 
on these populations is considered. A region in the ¢; — €2 plane where 
equilibrium points exist, is studied. .The stability properties of these equi- 
librium points are investigated. 


1. INDRODUCTION 


The subject of the effect of dispersal of populations isa topic of considerable 
ecological interest'-*. Holt* has considered a two patch model and a migrating predator 
from an optimal habitat selection point of view. Hasting® focused on spatial diffusion 
but also had a two patch model and showed the stabilizing effect of high dispersal rates. 


The view of the problem taken in this paper, as in Freedman and Waltman! and 
Freedman et al.?, is that “pressure” to disperse is given as a (monotone increasing) 
function of population size but that dispersal is inhibited by the difficulty of leaving the 
habitat which we, in turn, think of as surmounting a “‘barrier’’. It turns out that the 
more reasonable parameter is inverse barrier strength. This view appeared in Freedman 
and Waltman'. When this (vector) parameter is. zero, dispersal is impossible (the 
barrier is infinite) and each population grows to its carrying capacity. 


By the implicit function theorem, Freedman and Waltman' showed that, for 
small values of this parameter, the equilibrium was continuous and they approximated 
this equilibrium as an expansion in the parameter. This was done in Freedman and 
Waltman! for two habitats and a common barrier strength. 


Freedman et al.? discussed both the two habitats and n-habitats cases each per- 
mitted to have a different level of difficulty in its “escape” barrier. In addition, once 
the population has left its present habitat it may not successfully reach a new one (pre- 
dation harvesting, or for other reasons). In the analysis, by Freedman et al.’ they 
regarded the probability of a successful transition between habitats as given and ana- 
lyzed the question of the existence of the equilibrium as a function of the inverse barrier 
strength. Under reasonable biological hypotheses and one technical hypothesis, they 
determined the region exactly for two habitats and, in general case of n habitats. 
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2. THE MODEL AND EQUILIBRIA 
In this section we shall consider the case where a population is able to rer 
among 2-different habitats at some cost to the population in the sense that the 
probability of survival during a change of habitat may be less than one. This situa 
tion is described by a system of two prey and one predator of the form 


os eo a xi(1 Me | - By X,Y — € Xy + €2 Pas Xe 


sat cas a Xa(1 a* = | — Boxey — €2 X_ + €1 Pr2Xy 
1 


W—-—~+,-—— 


y =y(—y + 8, x1 + 82 x4) 
with Piz + Po; < 1 and x% (0) > 0, i = 1, 2 and y (0) > 0. 


Wa PSN 


x; represents the same population (prey) in the rwo different habitats; y is a pre- 
dator feeding indiscriminately on the two prey x1 and x2; B, and 8. measure the feed- 
ing rates of the predator on the two prey x, and x, respectively; is the death rate of 
the predator; 5, and 6, the conversion rates of prey to predator; €, and €, not necessarily 
small, but positive and represent inverse barrier strength in going out of the first habitat 


and the second habitat; and py, is the probability of successful transition from ith ha- 
bitat to /th habitat (where i + /). 


System (2.1) has been discussed by Freedman and Waltman! in the special case 
Piz = Px = | and ¢; = €, = € (the case of a common barrier strength) and for suf- 
ficiently small positive € (the case of strong barrier). In our paper ¢, and «, are positive 
but not necessarily small. We seek to find the region in the e, — «, plane where an 
interior equilibrium exists and to determine its stability properties. Also, the system 
(2.1) is discussed by Freedman et al.* in general case but without predator. 


Let E* (c;, €2) = (x¥ (a, &), x} (€1, €2), »* (€,, €)) denote an equilibrium point 
in positive octant, if it exists. Then, we have the following theorems: 
Theorem 2.1—E* (0, 0) exists if : 


cre We te. artes ee 
B, 3, ky 3 By Bs = B, 8; ky 
and (2.2) 
8:k, + Ske > AgeD 
where 
* (0,0) = FCO. ka ~ op 8 Bae + oa Bu) 7 
&, By Oy ky + ae By 8, k, ar: 
x* (0, 0) ——K2 (% By 31 ky &; Bo 3, k, + oy B, y) 
2 > <2 oy OL Soke ar eee eee 
%1 Bp 0, ky + a B, d, f 
y* (0,0) = 21% (8: ki + 32k, — y) 
01 B, 9p Ky + a By a ky j if 223) 
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Theorem 2.2—(a) E* (0, 2) exists if : 





B, 3; k, ay 3, ky } 
B» oF k, i Xo | % > | ! 
and ‘ 
ee ey. Feria .s bs ee 
2 B, a B, + Gy 5, ky | J ane 4) 
(b) E* (e;, 0) exists if : 
Bodske — af Sek, | } 
8, 5 ky Qy af 
and , 
ies = B sae re oY Z| 
: em B, 8, 3. k J . 28 Be) 
Theorem 2.3— E* (é,, €,) exists if either 
7 ; 
0<a<a(1—-F Joro<q<a(i—s?). 2.6) 


Proof of Thorem 2.1—Theorem 2.1 is stated in Freedman and Waltman! and its proof 
is a routine algebraic manipulation for solution of the algebraic system 





x 
an(1-~)-Bay=0 | 
| 
X2 
a X2 ( ks )- Bx» =0 
8, x; + 5, Xo = Y. J nee te & 


Proof of Theorem 2.2: Ate, = 0, equilibria are solutions of following system of 
equations : 


wy 
ami (1-2 ) —B; X,Y + & Po X2 = O 


a, (1 — 7) —bey —-@ = 0 


8; X1 + 9, X%. = ¥- 


From the second equation, we have 


ae et a 


“a(2:8) 


A, Xo 


Beets .(2.9) 


and from the third equation, we have 


= 62° %s 
5 = _t 3, fe . +++(2.10) 


Substituting x, and y in the first equation yields : 
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—a, Br ko (y — 8, x2) (81 ki — y + 8, x2) + Bi 8, kx (y — 8, x2) 


[(%. — €2) ky — d, X2] — fro k, ke ét Pos €2 X2 = O, 


this takes the form 
ax) + bx. +c¢=0 


where 
a = [a, By 6, ky + a, B, 8, ki] 8, 
— b = y [as B, 8, ky + a, B, 8, ky] + 32 ky [ag Bi 8, K, — a; By dy k, 
+ a B. ¥] + 8, k, ky €, [B, 8; Po: — B, 32] 
and 
C = [2% B,—a,B,] 51 ky key + a, By ky y? — Br kik 8, 7 ex 


From (2.10) for x1 > 0 must be x, < ¥/5,. 
Let f(x.) = ax} + bx, +c, then 


4, 2 
f ) Saas sb +e. 
5? %. 


Substituting and simplifying we get 


By 8 k, k, P2, &g 


f()-- ms lei: 
rar 


f (0) = c. Thus 0< x, < 1/85. 





And 


Ifc > 0, that is if 
a1 Y ’ 
sacteelae 


a a 
€5.<5 ieee Pak potest > 
‘ A ( Be By Be) 
From (2.9), we have 
Boke y = (a, — €,) ky — a x, 
> (% — &) k, e by (2.11) 
2 
> ok BY ag (oy my 
33 oot eee toe) 


a1 Be ks i Bekoy as y 
By By 33 ky 3, , 
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Thus, for »y > 0, must be 


ho 
ke a 
Yl ate t AG ERE 


or 





B, 3; k, Oy 5, k, 
Gtukor fx, ¥ Sa}. 


This completes the proof of (a), and the proof of (b) is analogus. 


Proof of Theorem 2.3 
This is the case where «, > 0 and «, > 0. Equilibria are solutions of the follow- 
ing system of equations: 
x, 
ax (1- = )-amPepan = Biny : 
1 
pe 4 
ay X2(1— 2 ) — exe tes pie x = Ba wey | 
2 
1 ant2z3) 
§,X%1 + 8,x%. = ¥. J 
Substituting the second and the third equations in the first equation of system (2.13) 
and simplifying we have 
2m | go (1 ba) — Me (y — 8, x,)? — str — 82%) 
0) 8? k 1 
1 


1 








€, B 
+e Pam |= ee Xo [ky — Xa] [ ¥— 8%] — 5 Xo 


Cea EL €; By Pig (y — 35 x,)*. 
5? 


1 


After simplification, we get the form 


ax, + bx, +cx,+ D=0 


where 
Gy B, dg ho By 5, 
, K, kK. 


b= a, B, 5; + By 8, €2 + By 5, Pio & _ B, 5, €; —ae B1 8; 


oie. « 
2a, Bev ae Boy Be O, Poy €2 
~ k, 3, ke 8, 
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2 
2 &, By dy ¥ al | GB, 8) y_ B. 3, Y «1 at, Bo ¥ 
saat ae = age a dy Es 8. 8, ky 
=a eee 2B, Y Pin 
2 
and 
ite Bi ¥* Pio & 
do ; 
Let 
f(x) Sax, + bxi + cx, + D. 
Then, we have : 
f(0) =D>0 
and 
3 2 
t(-—-) = at —— b+ ee 
= 53 5? do 


2 


Substituting and simplifying, we have 


¥ 8, oe 1? Poy €2 
t(4)=- ——_—_———_ < 0, 
53 


2 


Then, 0 < x, < y/6,. Similarly, 0 < x, < 7/9). 


From the second equation of system (2.13), we have: 
Bmy =| a —€ — | m+ Y Pie &1 a” 8, Pie €4 Xs 
k, 83 3; 





a * 
> E —- €&— 5k. : |x + a — do Bis €) 
= 1 1 2 





cad dae! eres lites st x 
d2 k, ae 


AX ¥ 
Then, B2y > a —«, — FE . Thus, for y > 0, must be 





2 < a,/ eee) 
dk, 


Similarly, — —!_. is 
Marly, €; > a, E 3, i ] - This completes the proof of Theorem 2.3. 
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3. STABILITY 


Having established the existence of an equilibrium E* (e,, €,), we proceed to ex- 
amine its stability properties. In fact we have the following theorems : 


Theorem 3.1— E* (0,0) is asymptotically stable. 

Theorem 3.2—(a) E* (0, €2) is asymptotically stable if ®: > 8.: 
(b) E* (e:, 0) is asymptotically stable lf B. > Bi. 

Theorem 3.3—E* (€:, €.) is asymptotically stable if 8: = £,. 


Proor : The first step is to compute the variational matrix of E* (€,, €2), which 
takes the form 











f x x ) 
| a (1 c =] =p y er. €7.— Bees €, Poy — Pa %1 | 
1 
V(a,a) = ; “e | 
ie «:) €) Pyo % (I— (x,/k,)) =e yy — ts — = — By X2 | 
2 
| 
t Sup dy 0 J 
or) GEN 
Making use of the system (2.13), the variational matrix takes the form 
( 9 
a re Sh | 
| 
V (€,, €&) = x h, Xe | 
(4 &) €: Pin — €1 Pie nr : 
| 
L ay Sey 0 J (3.2) 


The corresponding characteristic equation is given by 


x x ho xX x 
[> a a5 + €5 Poi | [a(a + gaa + €)Pro “| i BSny | 


2 


— €, Pro [€y Pr A — Bi 3, x; y) + 8, Y [eg Po, Bs X2 + B, x; 





x *. 9 
x (a+ ea +4 Py 2) | =o. (3.3) 


Equation (3.3) can be written in the form 


AS + a, A? + a,A + a3 = 0 





where 
a1 X1 he Xo x1 x2 
y= + + ey Py —— +2 Pie 
: k, k, X xy 
2 
€, &1 Pio x, 
Oy Ao Xi X ee 
a, = Se + Bi 8 ary + Bs Bs Hey + k, x2 
1 Ke 


(equation continued on p. 212) 
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€y % Pr x3 
ch ks X1 
and ne 
9 x 
a, = B, sy ae + €1 Pre 1) + BeBe x» k, 








+ €) Py2 42 | 4a B, 32 Pre x1 V + €2 B, 5; Poy Xo Y- 


Since both a; and a, are positive, then utilizing the Rough-Hurwize criteria, we get 
E* (e,, €,) is asymptotically stable if a,a, — a3 > 0. After simplification we have: 





3 
€, “1 Pie x? E> ho P21 Xs 
G1 ho xy X2 — 
ay a», a a) ; a 


k, ke k, Xo Bs ke x1 


2 


+ 3 aoe TEAR — Bo) €2 Poi X2 ] 





2 
= doy [ae “tr (Be ome B;) €) Pi2 x; } ee i te | 
2 
From (3.4), it is clear that if «, = ¢, = 0, then @, a, = a; > 0 and hence E* (0, 
0) is asymptotically stable. This proves Theorem sale 


Also a; a, — a3 > Oife, = Oand $: > 8: or¢, = O and 8, > B This proves 
Theorem 3.2. 


Finally, a, a, —a, > 0 if Bi = 8, for every positive values of €, and ¢, and this 
proves Theorem 3.3. 


Corollary —From the proofs of Theorems 2.3 & 3.3 it is easy to deduce the ulti- 
mate sizes of the two populations when E* (e,, ¢,) is asymptotically stable; for the prey 
x; (i = 1, 2) must be 0 < x; < y/3; where 5, x, + 5, x, = y that is x, x, belong 
to the open straight line segment AB where A (y/5;,,) and B(0,y/5.) Finally, for the 
predator y must be 





l ay vv oa Wie 
9 Se a max (a €&.- — , e— € — ~ 2 ) 
IS) vik, 35 k, 
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(Received 4 August 1987; after revision | August 1988) 


Itis known that most natural applications give rise to mechanisms of 
indirect actions, where the decisions in the control function 4 are shifted, 
twisted or combined before affecting the evolution. An example of this 
delayed action is the model defined by 


x(t) =e (x(t), ud) + f BU, 9) ult +5). 
—h 


We give sufficient conditions for the null controllability of such systems 
when the controls are limited in the sense that they lie in a closed unit ball 
with zero in its interior. These extend known results. 


1. INTRODUCTION 


Consider the nonlinear control system with distributed delays given by 


HD =AGXOUMXOF+LGEENUE+S) 


+ f(t, x (t), u (t)) ania) 


or 
£(t) = g(t: x(t), 4 (1) +f (d, B(t,s))u(t +s) (1.2) 


where x € E",u € IB =L ((a, 5], E™), the space of functions measurable and essen- 
oo 


tially bounded on finite intervals. It is assumed that the system (1.1) (= (1.2)) has a 
unique solution x = x (t) = x (f, fo, Xo, ¥), for each admissible control u € JB, and 
each initial state x (to) = Xo. 

The controllability of system (1.1) has been studied by Klamka‘”®, where he 


showed that under certain conditions on the matrix functions A (t, x, u), f (t, x, 4), 
B (t, s), system (1.1) is globally relatively controllable if 


inf det W (to, 43 Z, Vv) 2 ¢. 
Zu E Cam [to, tJ 
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Here 
W (Gy tii 2s v) 


is the controllability matrix (to be defined later). In this method, fixing the variable 
arguments of the matrix A and the function f by (z, v) © Cam [to, t:], he arrived at the 
linear dynamical system 


z(t) = A(t, z(t), v(t) x (0) +) (d, B(t, s))u(t + 5) 


+ f(t, z(t), v(t). 
He then applied Schauder’s fixed point theorem to this linear system to obtain condi- 


tions for global controllability of (1.1). 


In our own case we shall treat the equivalent system of (1.1) given by (1.2), 
where g (ft, x (t), u(t)) is nonlinear. Weshall employ the Alek seev-type variation 
of constant formula given by Khanh? to obtain an integral equation for system (I.2). 
We shall then use this integral equation to derive necessary and sufficient conditions 
for the null controllability with constraints for the linear variational system of (1.2). 
With this result, we shall then provide sufficient criteria for the null controllability of 
(1.2) with constraints. Indeed, we shall show that if 


(i) g(t, 0, 0) = 0; 


(ii) the system 
*(t) = A(t, x (t)) + B(t) u(t) + f(d, Bit, s))u(t +s) eet leat 


is controllable, 
where A (t) = D, g (t, 0, 0); B(t) = Dy g (t, 0, 0); 
(iii) the system 
* (t) = g(t, x (t), 0) : .. (1.4) 
is uniformly asymptotically stable, so that the solution of (1.4) satisfies 


LCE) Vissi] 9. H80=  eres 0, ¢ > 0, then (1.2) is null controllable with constraints 


Provided & Satisfies all smoothness conditions for the existence and uniqueness of 
solutions, 


We shall further apply our results specifically to the systems 


&(t) = g (x (t), u(t)) + Re Ciu (t — hy) salts) 


#(t) = Ax (t) + Bu(t) + na Ciu(t — hy) pith 
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where 
A = D, g (0, 0); B = D, g (0, 0). 1) 
Our results will show that if 
(i) g (0,0) =0 
SOR PNR eh A eens AM BR, ey oon cig At Cyl ee nt 
and 
(iii) the system 
#(t) = g(x(#),0) | .-.(1.8) 
is uniformly asymptotically stable then (1.5) is null controllable with constraints. 
Sebakhy and Bayounmi’ considered the system 
K=A(t)x(t) + B(thu(t) + C(t) u(t — A). 


This system is totally linear. Their results show that this system is null controllable if 
and only if 


rank I (ft t,) =" 


where [' (fo, t,) is the controllability Grammian. This result cannot be applied to 
systems of the form 


k(t)= g(t, x (t)) + B(t)u(t) + g(t) X(t) 
where g and fare nonlinear. Thus, our results complements and extends those of 
Sebakhy and Bayounmi’. 


Khanh? considered the nonlinear system 
z(t) = g(t, x,u(t)) + Bt, x(t) u(t) +f, x (0), u (t)). 
nonlinear though, but without control delays. Thus, our results extend those of 
Khanh? to include systems with delayed controls. However, our results are of a special 


nature in the sense that it considers controllability to the origin of E”, a results that 
is basic in the study of optimal control problems. Klamka** considered the following 


systems 
0 
& = A(t) x (t) a | (dH (t, s)}u(t + s) 
a linear system, and 


#() = A(x) x) + f GHGS) xu +s) 
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nonlinear. He gave sufficient conditions for controllability of the above systems when 
the controls are unlimited. Our result specialises the results of the above to control- 
lability to the origin of E" with limited controls i.e. with constrained controls. 
Yamatoto® and Dauer' gave results for the complete controllability of quasi-linear 
systems of the form 


& = A (x,u,t)X + B(x, u, thu +f (x, u, £). 


when the controls are unlimited. Our result strengthens and extends the above results 
to include systems with control delays. 


In section 2, we give the basic notations and define the system equations. 
Sections 3 gives the controllability and stability theorems for both the linear varia- 
tional system and the nonlinear system (1.2). 


Chapter 4 deals with the applications of the results to system (1.5). 


Throughout the paper, we exploit the properties of the ‘Reachable’ sets in deriv- 
ing our results. 


2. NOTATIONS AND PRELIMINARIES 


Let E denote the real line and J = [to, 4,], an interval in E. Fora positive 
integer n, we denote by E" the space of real n-tuples with the Euclidean norm denoted 
by |. |. If / is any interval of E, the usual Lebesgue space of functions, measurable 
and essentially bounded from J to E” will be denoted by L.. (J, E"). Mim will be used 
for the collection of all real n x m matrices with a suitable norm. Leth > 0 be given. 
For functions u : [t, — A, to] > E", t © [t, t,], we use “u; to denote the function on 
[—A, 0] defined by wu (s) = u (t + s), for s € [—A, 0). 

We shall consider the nonlinear system 


4 (1) = A (1, x (1), u(t) x (1) + Vd, B(t, s)) u(t + 5) 


+ f(t, x (t), u (t)) ad A) 


or 
4) =8 (x (0,U00) + FBS) u(t t 3) (2.2) 


satisfied almost everywhere on [to, t:], where the integral is in Lebesgue-Steiltjes sense 
res respect to s, x (rt) E E*, u € Leo ([to, th], E”), B (t, s)isannm X m matrix-valued 
unction absolutely continuous in S for each fixed t, and of bounded variation ins on 
[—h, 0] for each t € [fy, ¢,]. We shall assume that g (t, x (t) u (t)) (nonlinear in 
general) is continuous with respect to its arguments and continuously differentiable 


with respect tox. We shall als 
; © assume toat the functions i 
regarding the arguments. OS ee 
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Throughout the sequel, the controls of interest are JB = Lop ([to, f)], E™) 


IU © Leg ([h, th), E™), a closed and bounded subset of 1B with zero in its interior 
relative to 1B. 


Definition 2.1— The pair z; == {x (t), u:} is said to be the complete state of the 
system (2.1) or (2.2). 


Definition 2.2— System (2.1) is relatively controllable (or controllable) on 
[to, t,] = J, if for every complete state z; and every vector X,, x, in E" there exists a 
control u(t) defined on [fo, fi], such that the corresponding trajectory of system (2.1) 
satisfies the condition x(t:) = x,, with x (t) = X 9. Following Sebakhy and 
Bayounmi’, we give the following definition of null controllability. 


Definition 2.3— System (2.1) is said to be null controllable at tf = f1, if for any 
initial state {xo, U,} on [t, — A, to], there’ exists and admissible control u(t) € /B, 
defined on [f, t; — A] such that the response of the system is brought so the origin of 
E”" att = tf, using the control effort 


u(t) = {u (t) on, [fo, t1 — Al} 
= 0 on, [h — A, 4). 
That is if given the initial state {xo, u,} on [to — A, fo] there exists a control u (t) € IB 
defined on [f,, t; — A] such that x (t,) = 0. 


System (2.1) is null controllable with constraints at tf = ¢, if for any initial state 
Xo, Ug ON [t) — h, to], there exists an admissible control u (t)E IU, defined on 


[t,, t; — A] such that the response x (f) satisfies x (ti) = 0, using the control 


{* (t) € 1U on [to, t; — a 
AUS Bor 0 € JU, on [t, — fh, th). 
Let the system 

x(t) = g (t, x(t), u (1) nah 219) 


have a unique solution x = G = G (1, fo, u, 1), for each admissible control u € /B 
and each initial state G (fo, fo, ¥, u) = »y E E*. If we now set 


Ft, lo, Ys u) = Gy (t, to, J, u) ...(2.4) 
and 
K (t, x, u) = By (t, x, u) Ent oy 


then it is not difficult to see that the following relations hold for the n X n matrix 
functions 


Fi (t, s, y, u) = K(t, G (t, 5, yu), U(O) F(t, 8, Ys Ws F (los fos uy Y) = Inn 


218 JERRY U. ONWUATU 


F (t, lo» u, A (to)) = t ee j K (s, G (s, los uy (to), u (s)) 


F (5s, to, V (to), u (s)) ds. alos0) 


Thus, Alekseev-type variation of parameter formula for (2.2) (= (2.1)) is given via the 
result of Khanh? by 


x (t, bos Xo, u) a G (t, bos Xo» u) + f FAr, 5S, xX (s), u (s)) 


0 


[y (d, B(t, s)u(t + s)) ds}. so(2:%) 


When we consider the solution x (t) of (2.2) given in (2.7), we note that the values of 
the control u (t) for t € [t) — A, to] enter into the definition of the complete state 
z(t )). Thus, the last term of (2.7) must be transformed to take care of this by inter- 
changing the order of the integration. Now, using the unsymmetric Fubini theorem 


(Klamka‘), we have the following : 
t 


o 0 
X (t, lo, Xo, U) = G (t, bo, Xp, u) + J dB, (jf F(t, I —s,x(l —s), 
-h tts 
0 
0 
u(l— s)) B(l — 5, 8) us dl) + | dp 
-h s 


CV F@I~ sx 9),u(l— 9) BU — 5,3) uO dl 
(2.8) 


where the symbol dg, denotes that the integration is in the Lebesque-Steiltges sense 
with respect to the variable s in the function B (t, s). If we denote by 


Bil, s)forl St,sEE 
} «aoe 


B; (i, s) = | 
0 forl>tsEE 


then (2.8) yields 


0 9 
t, lo, ’ —~ 
X(t, to, Xo, u) = G (ft, to, Xo, u) + Le dp. ms F (t, l—s, x (l—s), u(I—s)). 


B(l— s,s) dl) + J dof F(ty1 — 5, x,u) B(l—s,s)u() dd) 


...(2.10) 
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t 
0 0 
X(t, fo, Xo, u) = G (ft, t; xXo,u) + J ds,( J F(t,] —s, x, u) 
—h (rs 
0 


(B(I — 5, 5) ur, dl) + f F(t, 1 — s, x, u) 
t =h 


0 
d, B; (1 — s, s)) u (1) dl. REY 
We now define the n X n controllability matrix of (2.2) (= (2.1)) by 


W (to, fi) = J 2 F (t;, 1 — s, x, u) d, B, (J — s, s)) 


( i F (t,,1 — s,x,u)d, B (I — s,s)? dl. 2.123 
= 1 


Where the symbol T denotes the matrix transpose. Note that W(t, t;) is symmetric 
and non-negative-definite. 


Definition 2.4— The reachable set of (2.2) at time ¢, using L.,.control is the 
subset of E” given by 


IP (ty, to) = [f F(t, s, x (s), u (s)) J (d, B(t,s)u(t+ s)) ds: u © JB}. 


0 
In a similar manner, we define the constraint reachable set of (2.2) as 


IR (t,, to.) = [ TFC, 5, X (s), u (s)) i (d, B(t, s)u(t + s)ds: u € IU). 


0 


3. CONTROLLABILITY THEOREMS 


Consider the system 


x(t) = A(t, x (t), u(t)) x (t) + I (a. Bf s)) uit +3) set stl) 
or 
4 (t) = (t,x (1), u()) + i (d, Bit, s))u(t +). 8.2) 


Lemma 3.1— For each fixed (t), 50, Xo, uo) © E* xX E* x E"® X Lo, let 
Dx (to, So, Xo, Uo) denote the partial derivatives of x with respect to its last two argu- 
ments at the point (fo, 50, Xo, 40). Then, for every (h,u) € E" x L.., we have 


D (to, So, XO, Ug) (A, u) i A (to, So, X05 Uas h, u), 
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where the mapping t > A(f, So, Xo, Mo, 4, u) of Einto E” is the unique absolutely 
continuous solution of the linear differential equation 


z (t) = D, g (t, X (t, So, Xos Uo), Uo (£)) Z (t) 


+ Dy g(t, X (ty So» Xos Mo), Mo) V(t) + I (d, B(t, s)) v(t +3). 


A ee WEY 
Also, for every W € L., we have 


Dd, x (to, So, Xo, Uo) ve v (to, So, Xo, Uo, w) 


where the mapping f — Uv (f, 59, Xo, Uo, W) is the unique absolutely continuous solution 
of (3.3) satisfying the initial condition (So, 0) i.e. z (59) = 0. 


Proor : Let D, be the partial derivative of x (t, 59, Xo, u) with respect to u. 


Then, 


D,, X (So, So, Xo, Uw) = 0 in [—A, 0] 


‘ 

D, x (t, So, Xo» u) pm J D, g (s, x (s, So, Xo, u), u (s)). D, x (s, So, Xo, u) ds 
5 
0 


“+ j Ds g (8, X (8, So, Xo, U), u(s)) ds 


+1 Fd, Bs) ds 3.4) 


Differentiating with respect to f, we obtain 
d 
ar Du X (ts So, Xo, ¥) ¥ = Dy g (Ss, x (5, $0, Xo, u), u). Dx (8, Soy Xo, U) ¥ 
0 
+ Ds g (Ss, ¥ (8, $5, Xo, U), BE | (d, B(t, s)) v. 
-h 


Substituting the limits of integration in (3.4), we obtain «»(3-5) 


18 ye x (So, So, Xo, u) i TPs x5 v= 0 
proving the second assertion. To prove the first part, we note that 

DxAt, Sox tenth) kb) ae Ds X (t, So, Xo, U) + Dy x (t, So, Xo, u) v. 
But, 


t 
Ds X (t, So, Xo, u) = T+ J Dy g (5, x (5, 50, Xo, u), u (S)) 
0 


X Ds x (5, 59, Xo, u), u (s)) ds 
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D, X (8, &, Xo, 4) = f Ds g (5S, X (S, So, Xo, u), u (s)) ds 


+5 UF Gi, BU, Mas 
Therefore, 


Dx (t, So, Xo, u) (A, v) =h + f D, g (s, x (s), u(s)). Ds x (S, So, Xo, u) hds 


0 


‘ 
+ J Dz g (s, X (5S, 5, Xo, u), u (s)) vds 
‘9 


+f if (d, B(t, s)) vds. (a6) 


Taking the ¢ derivative, we obtain : 
d 
“dt Dx (t, So, Xo» u) (h, v) = 0 =e D, g (s, x (s, So, Xo» u), u (s)) 
xX Ds x (S))h + Dag (Ss, X (S, So, Xo, U), u (s)) v 


oe i (d, B (t, s)) v. 


Moreover, from (3.6), we have 
Dx (50, 50, Xo, u) (A, v) =h+0+4+0+40. 
Hence, the Proposition. 


Definition 3.1— Let X, Y be Banach Manifolds of class C' and leth: X¥ > Y be 
a C! mapping. A point y, in A (X) is called a normal value of / if there exists at least 
one Xo € A-! (yo) such that the differential dhx, is a split-surjective linear mapping. 


The next Lemma gives a necessary and sufficient condition for the derivative 
D,, x (t, So, Xo, Uo) to be a surjective linear mapping, in which case the point 
x (t, So, Xo, Uo) is a normal value. 


Lemma 3.2— Given %) € E",u € Lo. For t, > 8 +h, lett € [9 - fi]. Let 
u —> X (ty, So Xo, U) be the mapping T: Ln. > E£", defined by 


T (u) =a (4, Sos Xo» u) 


d n 
where x (f1, So, Xo, ¥) is a solution of (3.2). Then DT (u) = oi T (u): Loo > E" hasa 


continuous local right inverse (is a surjective linear mapping) if and only if the linear 
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variational control system (3.3) of (3.2) along the response f — x (f,, So, Xo, u) is con- 
trollable on [s9, fi], t; > 59 + A. 


Proor: For t € [5 9, t;], t; > So + Ayu E Loo, lett > z(t, to, Xo, u) denote the 
solution of (3.3). Clearly, the linear variational control system of (3.2) along the 
solution tf > x (f, 59, Xo, u) is controllable on [sy, #,] if and only if the mapping 
u —> Z(t, So, Xo, 4) is surjective, f E [so, f,]. But by Lemma 3.1, 


Z(t), So, Xo, U, vy) = D, x (h, Sos Xo» u) = DT (u) y. 


Therefore, the marpirg v > z (t,, $0, Xo, u, ¥) is surjective if and only if DT (u) is 
surjective. The following proposition will te mzede use of in the next theorem. 


Proposition 3.1 (Surjective mapping theorem) (Graves)— Let U be openina 
Banach space X. Letf: U-~ Y be aC’ mapintoa Banach space Y. Let x) € X. If 
f' (Xo) is surjective, then fis locally open in a neighbourhood of Xo. More precisely, 
there exists an open neighbourhood V of x, contained in U having the following 
Lang® property : 


For each x € V, and open ball Bx centered at x, contained in V, the image /( Bx) 
contains an open neighbourhood of f(x). 


Remarks 3.1: This is another form of the implicit function theorem given by 
(p. 193), 


Theorem 3.1— The system (3.2) is controllable on [so, t], t > 59 + A, whenever 
the system 


2 (1) = Deg (1,0, 0)z (1) + Dy (1,0, 0) 9) + f (d, BUt, 5)) v(0-+5) ds 


...(3.7) 
is controllable on [s9, t], t > 5s) + A. 


Proor : For system (3.2), let vu € L.., and x (so, Xo, u) be its response. Let 
T: Loo E” be the map defined by Tu = x (t, 50, 0, u). 


It follows from the conditions on g and (2.11) that T (JU) = IR (1, 59). Suppose 
system (3.7) is controllable, then by Lemmas 3.1 and a Bit J 


d 
DT (0) = du [T (u)Juxo => D, », (1, So; 0, u) | 4u= 0 


is a surjective linear mapping of L,,—> E”, 
open. Hence, there is an open ball Pig Gt & 
containing zero such that B, C 2 (B>): 


Therefore, by Lang® (p. 193), T is locally 
oo Containing zero and an open ball B, CES 


Since /U contains an o ini 
pen ball containing zer ; : 
such that B, C T (Bp A IU), Beene te Saar 
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Therefore, 
so that 0 € Int JR (f, 59). 


Remarks 3.1 : The above theorem assures us that if the linear variational con- 
trol system of (3.2) corresponding to the zero solution is controllable, then 0 € Int 
TRU ft). 


Definition 3.2— The domain D of null controllability of (3.2) is the set of all 
initial points x, € £", for which the solution of (3.2) with x (t)) = Xo, satisfies 
x (t;) = 0 € E*, at some f¢, using U € JU. 


Theorem 3.2— In system (3.2), assume 
(i) g(t, 0,0) = 0, t > So 
(ii) system (3.7) is controllable on [s9, t],f > so +h 
and 
(iii) the smoothness assumptions hold for g 
then 


the domain of null controllability of (3.2) contains zero in its interior. 


Proor: The variation of parameter formula for system (3.2) is given from (2.7) 
as 


ea) = O Ue taxa Ww) + J F (t, s, x (s), u(s)) 


i B(t, s))u(t + s) ds. tae) 


We note that if there exists a u © /U such that the solution of (3.2) satisfies 
Mi facta Xe Hy = Xoy (Li, bo, Xo, 4) = 0, 


then from (3.8) we have 


0=G(t, tr %u) +f F(t,s,u, x (9) a, B(t,s)) u(t +s) ds. 
hy - 


By (iii), (i) 0 € D, the domain of null controllability of (3.2). Assume that 0 & Int D. 
Then there exists a sequence {x)}“i-, C E" such that xi > 0 asi > co and no x 1s 
in D, hence, x; ~ 0 for any i. Since x; & D x (1, fo, 1, u) ~ 0, for any t > % + h 
and any u € IU. Hence, from (3.8), 


~G(t, te xy) SHS F(t, 5,4, x(9))[ inG? B(t, s,)) u(t + s)] ds 


0 
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for any i, any tf > tf) +h and u € JU. Therefore, ¢; is not contained in JR (t, t,) for 
any t > tf, + h. Hence, we have a countable sequence {%), C E”, such that 

G1 > 0 asi co 

GC & R(t, to), for anyt >t +h 

C ~ 0 for any i 
hence 0 & Int /R (rt, t.), for anyt > t) + fh, 


But by (ii), the linear variational system of (3.2) about the zero solution is con- 
trollable, hence, by Theorem 3.1, 0 € Int JR (t, to) t > fy +h which is a contradic- 
tion. This contradiction shows that 0 € Int D. 


Remarks 3.3 : The above theorem says that, if the linear variational system is 
controllable, then (3.2) is locally null controllable with constraints. 


Theorem 3.3— For system (3.2), assume 


(i) g satisfies all smoothness conditions for the existence and uniqueness of 
solutions; 


(ii) g (¢, 0, 0) = 0; 
(iii) system (3.7) is controllable on (fo, i), t) > to +A; 
(iv) the system 
% (t) = g(t, x (1), 0) --.(3.9) 
is uniformly asymptotically stable so that every solution of (3.9) satisfies 
Il X(t, to, Xo, 0) || Sk | ene) ps tk >0,a>0 ---(3.9a) 
then 
(3.2) is null controllable with constraints. 


PROOF : By (i), (ii), (iii) the domain D o: null controllability of (3.2) contains 
zero in its interior. Therefore, there exists a ball B, € Int’ D. 


By (iv) every solution of (3.9) which is a solution of (3.2) (with u = 0) satisfies 
x (t, to, Xo, 0) > Oast > cc, hence, at some ft, < eo, we have x (ty, fo, Xo, 0) = Xo 
€ #, C Int D. Therefore, with ,X_ as initial data there exists u & JU and some 
!, > t,, such that the solution x (t) satisfies x (t,, ‘1, X2, u) = 0; proving the theorem. 


We now give necessary 


and sufficient conditions for the controllability of the 
linear variational system (3.7) 


Now, consider the linear system 


#() =A) x(t) + But) a (d, H(t, s))u(t +5), (3.10) 
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Following Klamka‘, the unique solution of (3.10) is given by 


t 0 
X(t) = F(t, to) x (to) + J F(t, 7) [ J (d, H(t, 5)) u(t + 5) + B(t)] de 
ty = 
pee OA 
where F (t, s) is the fundamental matrix for the homogeneous system 


*(t) =A (t) x (2). meds i2) 
Following Klamka‘ and using the unsymmetric Fubini theorem, we can transform eqn. 
(3.11) to the following convenient form : 


t 


eye Rit) x (i) 4- Fda | joFt ane Oe ee WRENN 


A [ j F(,1— s) dH (I — s,s) + F(t,]) B()u(I) dl}. 
t -h 


0 


..(3.13) 


We now define the n x n controllability matrix of (3.10) as 
ae = 0 
W, (to, t,) = J tol F(t, t—s) dH (i= 5; MU F (ti, t—s) dH (t—s, s)]? dt 
t —h =. 


0 
fs 
+ ) (F(t, 1) BOLE (4, 1) Bt) dt. ...(3.14) 
t 
0 
The proof of the next theorem follows exactly the same proof given in Theorem | of 
Klamka® and therefore will be omitted. 


Proposition 3.1— System (3.10) is controllable on [fo, ¢,] if and only if 
rank W, (fo, t,) = m3 t; > to. Aas) 


Remark 3.4: Here, we extend and complement this result of Klamka*® to system 
(1.1) in the following manner. 
The next theorem is a consequence of Theorem 3.3 and Proposition 3.1. Let 


k(t) = A(t) x(t) + B(thu(t) + Sd B(t, s))u(t + s) ...(3.16) 


where 
A (t) = Dz g (t, 0, 0), B(t) = Dag (t, 0, 0) (3.12) 


and let 


W, (to, ty) =sIS F (th, t — s) dB(t — s, s)] Uf F (t, t—s) dB (t—s, s)} 


0 
(equation continued on p. 226) 
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‘ 


+ [F(t ) BODIE (4, OF dt (3.18) 


where F (t, s) is the fundamental matrix for the homogeneous system 
& (t) = A (t) x (ft). ...(3.18) 
Theorem 3.4— For system (3.2), assume 


(i) g satisfies all smoothness conditions for the existence and uniqueness of 
solutions; 


(ii) g (t, 0, 0) = 0; 

(iii) rank W, (to, t:) = m3 ty > ty 
where W; (fo, ¢,) is given by (3.14); 

and 
(iv) the system 
%(t) = g(t, x (t), 0) (3.19) 

is uniformly asymptotically stable then 
(3.2) is null controllable with constraints. 


Proor : Immediate from Theorem 3.3 and Proposition 3,1. 


4. APPLICATIONS 


If we now specialize to the constant system with multiple delays in the control 
defined by 


*(1) = g (x(t), u(t)) + z Ci u(t — hy) (4.1) 


then the following results follow : 


it is known! that the system 


k 
t(t)= A + 
x (t) x(t) + Bu (t) 4 a Cru (t — hj) »+.(4.2) 
is controllable on [0, ¢,], 4; > max Vig..wers » Am} if and only if 
rank [B, ..., A") B, Cy, |... 4n- Cissy Ck, ..A9-) Cy] = 1M ... (4.3) 


, Remarks 4.1 : The above results of Sebakhy and Bayoummi? cannot be applied 
directly to the nonlinear system (4.1), without some modifications as will be shown by 
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Proposition 4.2. The linear variational system of (4.1) corresponding to the zero solu- 


tion is given by 


z(t) =D, g (0,0) z(t) + Ds (0, 0) u(t) +3 Bee hose ath) 


If we set 
A = D,g (0, 0); B = Dy g (0, 0) ..(4.5) 


and suppose that A, B are autonomous, then the next Proposition follows from the 
reasoning above. 


Proposition 4.1— The autonomous linear system (4.4) corresponding to the zero 
solution of (4.1) is controllable if and only if 


Pattee eee oy CE AY Cy) = 
where A and B are given by (4.5). 


The next Proposition follows from Theorem 3.4 and Proposition 4.1. 
Proposition 4.2— For system (4.1), assume 
(i) (4.2) is autonomous; 


(ii) g satisfies all smoothness conditions for the existence and uniqueness of 
solutions; 


(iii) g (0, 0) = 9; 
(iu) rank [Bj 24 Ae Bs Caron A*-1 C,] = n, 
where A and B are defined by (4.5); and 
(v) the system 
* (t) = g (x (0), 0) ...(4.6) 
is uniformly assymptotically stable then 
(4.1) is null controllable with constraints. 


Proor : Immediate. 
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Some structural considerations are further made on the previously intro- 
duced Finslerian metrical structure g x (x, y) = YAx (x) + Ayn (x, y). In 
particular, some interesting features underlying the field equations for 
Finslerian gravitational field based on this metric are investigated. 


1. INTRODUCTION 


As has been mentioned ina previous paper’, there appear two fields around 
each point x of the Finslerian gravitational field : One is the external (x)-field spanned 
by points {x}, which is nothing else than the Einstein’s gravitational field?, and the 
other is the internal (y)-field spanned by vectors {y}, which is regarded as the so-called 
internal space attached to x. The former is dominated by the four-dimensional 
Riemann (R,) metric yyx« (x) (x, A = 1, 2, 3, 4), while the latter is assumed, in general, 
to be governed by the four-dimensional Riemann (Ry) metric /;, (y) (i, 7 = 1, 2, 3, 4). 


From the vector bundle-like standpoint‘, the (y)-field may be regarded asa 
fibre at the point x of the base (x)-field and the total space of this vector bundle may be 
considered a unified field between the (x) and (y)-fields. The Finslerian gravitational 
field, therefore, may be likened to this unified field, presenting an aspect of eight- 
dimensional Riemannian structure (R;) dominated by the R,-metric Gag (X4) 
(X4 = (x*, y*); A = (x, i) = 1, 2, 3, ..., 8) (Miron’ and Miron and Anastasiei*). This 
R, structure has been reduced, in the previous paper’, to a four-dimensional Finslerian 
structure (F,) based on the Finsler metric gyx (x, y) (2.9) by means of a dimension 
reduction-process. These situations will be considered in more detail in the following. 


2. ON THE FINSLERIAN STRUCTURE 


In the previous paper (see section 3 of Ikeda’), we have taken account of 
inherent law internal vector y (= y‘) in the form 


p= x; (x) y/ my te Y 


where K; represents the rotation matrix. From (2.1), one typical intrinsic paralle- 


lism (i €., connection) of y such as 
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Syt = dy' + Ki, yt dx* (dy + ML dx*) = 0 (2.2) 


eK, 
Tate and dy! = 7‘ — K; (0) »’. The quantity N,. 
in (2.2) plays the role of nonlinear connection in the theory of Finsler spaces*® and 


also serves, in our case’, aS the mapping operator of the internal (y)-field on the 
external (x)-field (see below). 


has been proposed, where Ke = 


In our Finslerian field, which is regarded as the unified field mentioned above, 
the connection relations are prescribed by*"* [see (3.4) of Ikeda’] 


DV® = dV*¥ + FY Vadxt + Ch, VA dyk 


DV‘ = dV‘ + Fi. V4 dx + Ci, V4 dyt om ey 


where 


K i 1 
(Free Caen) 


i © 


= (29) 


denote the coefficients of connection. From (2.3), the following covariant derivatives 
are introduced : 


A K « ovr sy 
Ue = 3e + Fi.V4, Vie = aye Ee Tid Se 


V 


SVs i 
Vion = — + ie V), ad ~—_ OV 


(2.5) 
35x — aye ae Cr Px 


——-+--— 


Therefore, the base and dual base (i.e, the adapted frames) are set by 


iO (Beall Oe iba 
oC4 ~~ \Sxe ax® * ayi? ay) . | 


ve 


deA = (dx¥, By! = dy + NE dat) | 8) 
J 
From (2.6), the so-called decomposition factors are defined by 


4 =), 0), 4, = (8, — N\), 


mas lek 
B= (Ni, 8!), B4 = (0,3!) (2.7) 


~-—-~-—~ 


by which the unified metric of the unified field Gas is decomposed to, under the 
assumption that Gax = yak (x), Gi; = hi, (y) and Gyi = Gin = 0 
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Bx = AL AP Gas = yax(x) + Ni N! hy) | 
| 

Zu = At B* Gap = — Ni hy (¥) 4 ...(2.8) 
| 

si; = B By Gus = hi, ()). | 


At this stage, as has been considered in Ikeda’, if the ())-field, is compactified, then 
only the F,-metric gyx (x, y) (2.8); appears above the surface. This compactification 
is likened geometrically to the mapping process of the (y)-field on the (x)-field, where 


the nonlinear connection Ny plays the role of mapping operator. (In the previous 
paper’, one new operator ey (x), instead of Ny (x, ¥), has been introduced as the 


mapping operator). By this mapping process, those quantities such as $y', oe at 


Cr etc. are brought to the external (x)-field and the unified field presents an aspect 


of F,-structure based on 


Ban (X,Y) = yan (xX) + Aye (x, ¥); 
Ayr (x,y) = Ny Ni hi (y). 2:(29) 


In the following, we shall pay our attention to some structural features under- 
lying the field equations for this Finslerian gravitational field from the vector bundle- 
like standpoint based on (2.4). 


3. ON THE FIELD EQUATIONS 


Now, in our theory, the (y)-field itself is dominated by R,-metric Aj, (y), as 
mentioned in (2.8), so that its field equation may be written as 


M1; (y) — & M hi; (y) = 11; (YY) wi(G2) 


where M;, denotes the Ricci-tensor derived from the Riemannian curvature tensor of 
the (y)-field, M the scalar curvature and »,;, means the energy-momentum tensor for 
this case. Then, (3.1) must be mapped on the external (x)-field governed by R,- 
metric y,x (x) by means of the mapping process mentioned above. As the result, 
therefore, the field equation of the unified is given by 


Ry« — 3 R grx == T)K AT te fee} 


° K 
where Rx is the Ricci-tensor derived from the third curvature tensor Ry,» formed 


with F, ~ 


case, since gyx is given by (2.9), Ry« has such a special form as 


R = Rxj g** and t,x denotes the energy-momentum tensor. In our, 
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etl Ees: 
Ryx (F) = Ky« ({}) + Max (A) (3.3) 
where K x represents the Riemannian Ricci-tensor, which is formed with the Christo- 


i c= 0 i.e., 
to the definition M)« = Ni N! Mi,). In (3.3), even if Ry 0, Myr S 
M,,~0). This case corresponds to the compactification of the (y)-field considered 
by Gasperini®. 


Next, we shall take up the following field equation for the empty space : 


Sones 6 (3.4) 


K 
where Sy, means the Ricci-tensor derived from the first curvature tensor Syae 


defined 


by ee (cf. Miron and Anastasiei* and Matsumoto®). Concerning this, it has been 


shown’ that if Sy, = 0, then pee = 0 holds good, so that the field itself presents an 


almost Riemannian aspect, due to Brickell’s theorem’. That is to say, the Finslerian 
field with Sy. = Q becomes almost Riemannian. Therefore, (3.4) seems somewhat 
unsuitable from our Finslerian viewpoint. Further, following the definition of 


ARE the Ricci-tensor Sy, (= Syax) in our case based on (2.9) is calculated as 





Sv, = $ yk Eva. hva - oo 
rae ae oy™ dyk oy* oy” 


sac Sua) 


at least in the first order approximation with respect to Ayx. In (3.5), if Aye 
22 : : 
=} Sree. then Sy, = 0 (3.4) is identically satisfied. This is also unsuitable 


from a physical point of view®. 


The field equations (3.2) and (3.4) for the unified field may be reconsidered from 
the vector bundle-like standpoint as follows: In the total space, the adapted frame 
(2.6) has been set and the connection with the coefficients (2.4) has been introduced. 
This connection is made metrical for the metrical structure of the total space 


G = Gap dtA dt® = gyx (x, y) dak dxd 4 81; (x, y) Syt 


dy! acta) 
In this case, 


the compactification Process with respect to (2.8) is not taken into 


account and the frame (2.6) is assumed ‘to be suitably adapted to the conditions 


fyi = 81, = 0. By straightforward calculations, the following six kinds of curvature 


tensors in the total space are obtained through the Ricci-identities with respect to the 
convariant derivatives (2.5) : 
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K 


eee K i i 
=A BCD... Uses Ry yps P eg? Sy ’ ‘ie ). Se ee h 


jk? 


As to the Ricci-tensors, they are given by, from (3.7), 


K 


A eee ee, se 
AaB (= Rage) = (Rav = Ris Pi, =P 


Si; = S\,). ...(3.8) 


At this stage, following Mircon‘, we shall define the field equation for the total space 
in the form 


Rag — 4D Gap = TAs tae) 


where the total scalar 7% (= Gaz G4) is given by 7 = R+S (R = Ryx g*4, S = Si, 
g‘’) and t,s means the eight-dimensional energy-momentum tensor depending on 


(x*, yf). (3.9) is decomposed into the following four kinds of field equations* by 
(3.6) and (3.8): 


Ry« — $ (R + S) gx« Stk 


Py = 71, 


P*y; = T37 ...(3.10) 
Si, =a 4 (S 1 R) 8ij —— Thy. 


These equations can be obtained quite systematically without any special conditions 
on the spatial structure itself. Of course, they depend on the choice of adapted frame, 
so that different decompositions of (3.9) result in different equations. It turns out, 
therefore, that (3.2) and (3.4) are regarded as special cases of (3.10). 


In (3.10), the most essential feature is the mixing of the Ryx-component repre- 
senting the x-dependence and the S;;-component representing the y-dependence. This 
means the interaction or coupling between the miro- and macro-degrees of freedom. 
The role of S in (3.10), is compared to that of cosmological term, because S is consi- 
dered an eight-dimensional effect (cf. Kerner’), It should be remarked that in the 
ordinary generalized Kaluza-Klein theory®”°, such components as P?;, and P*y; (i.€., 1, 
and t;) are not taken into account from the beginning. It may be considered that 
(3.10) contains some more instructive information, so that its physical aspects should 
be investigated in future. 
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In this paper, we give sufficient conditions for the existence of periodic 
solutions of the nonautonomous equation (1.5). 


1. INTRODUCTION 


Consider the fourth-order constant coefficient differential equation : 


+ax+a,xt+ax+ax=0. (Ly) 
As it was shown in Ezeilo® that if 

a,>}a? V1.2) 
or 


gas <— 0; a, 5-0 ... (1.3) 


then the auxiliary equation corresponding to (1.1) has no purely imaginary roots what- 
ever. By the general theory; this, in turn, implies first of all that (1.1) has no periodic 
solution whatever other than x = 0, and secondly that the perturbed equation 


x4) + a, x + a; a Or eS a,x = pit) ...(1.4) 


in which p (~ 0) is any continuous w-periodic function of f, has an \-periodic solution 
subject to (1.2) or (1.3). 


In the literature, there are some extensions of (1.2) or (1.3) in one form or other 
of the existence result for eqns. (1.4) where some of a,,... @, are not constants’”’. 
The object of the present paper is to extend the result (1.2) for (1.4) to equation in 
which a, a., a3 and a, are not all constants. 


We shall be concerned with the equation 


4 f, Ox thO eth +h =pl,x% xx) .(5) 


Where fj, fo, fz, p are continuous functions depending only on the arguments shown 
and p is also assumed to be w-periodicin?. That isp (1, x, ¥, Z, u) = p(t+w, x, 
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y, z, u) for some w > O and for arbitrary 1, x, ), z, uv. We shall however require here 


that /{, (x) exists and is continuous for all x. 


We shall establish here the following theorem 7 
Theorem—Suppose that 
(i) there exists a constant a, > 0 such that 
| fe (v) | © @ for all y (1.6) 
a, = inf f* (x) > 4.43 at EEA 
x 


(ii) there are constants Ay > 0, Ai = O such that 


| p(t,x,y,z,u)| < Ag t+ Ai(l|y|+])2]1) soeld Op 
for all ft, x, y, z and wu. 


Then there exists a constant €) > O such that (1.5) has at least one w-periodic 
solution for all arbitrary f, and fs if A, < é>. 


It should be noticed that we get Theorem 3 given in Tejumola’ under weaker 
conditions if we take f/f: (%) = a, and fy (x) = a, x in eqn. (1.5). 


2. PRELIMINARIES 


Asin Ezeilo*, the proof will be bythe Leray-Schauder technique, with the 
equation (1.5) embedded in the parameter-dependent equation : 


x + phi (x) x +{1 — wa, + af x + ef) +(1 — pax 


+ fi) =bp(t,x,x %x),0<p<l. (21) 


Note that when » = | (2.1) reduces to the original eqn. (1.5). Also when » = 0 
it reduces to the linear equation 


- 


x) + ax +a,x=0 


which, in view of the condition (1.2), has no non-trivial w-periodic solution. Thus 
Theorem will follow from the usual fixed point considerations in Ressig er al.? if it can 


be shown that there is a constant D whose magnitude is independent of » (0 < » < 1) 
such that, if x (1) is any w-perodic solution of (2.1), then 


IXOI<DIXOI<D| XO1<D[xWl<SD  ...22) 


for all t € [0, w]. Note that the f-ran 


ge here may be replaced by [7, T + w] (arbitr- 
ary T) since we are dealing with a w 


periodic x (f). 


Before proceeding to the actual verification of (2.2) we shall introduce some 
notations. Throughout what follows, D’s with or without subscripts denote finite posi- 
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tive constants whose magnitudes depend ona, @;, Ao, f; and f;. The D’s are all 
independent of ». Finally aD without a subscript is not necessarily the same 


each time it occurs, but the numbered D’s: D,, Dy, ... retain a fiixed identity through- 
out. 


3. Proor oF THEOREM 


We shall take (2.1) in the more compact from 


MO pfs (2) x + few CX 4+ whe) thaw bh () = pp ee D 


(O< p< 1) peal) 
by setting 


Soyp (4) = (1 — #) Ge + bfy (4) 
Sap (X) = (1 — pw) aa (x + Bf, (2). 


In what follows in the rest of this paper x (¢) is an arbitrary w-periodic solution 
of (3.1). It will now be shown that x (t) satisfies (2.2) if Ai is sufficiently small. 


Our main tool in the verification of (2.2) for x (t) is the function, 


x : x x 


Piet lant: d: +2 ‘| fam (V) dv + XX +X fap (X) + | 
fa (y) ay. 
An elementary differentiation gives 
V(t) = Us + px p(t, x, x % x) ...(3.2) 
where 
x 


Uy ax + ¥ | fn OD dd + Soy 
0 


> x —a,| x| |x| + ax? 8.3) 
by (1.6) and (1.7). 
It is now not difficult to verify that 
U, > Dy (x? + x*) (3.4) 
for sufficiently small Dy. For, suppose for example that Dy < 1, then, by (3.3), 


U, — Dy ( (x? + x 2) > > (1 — Dy) x? — a,|x| | x1 + G@ — Do) x" 
> 4(1 — Dy)? YU, # 
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where 
U, = (4a, — a3) — 4D, (1 + a4) + 4D8. 


However, by (1.7), 4a, — a3 > 0 and so U, is strictly positive if, say 


0 < Dy < 3 (4a, — a3) (1 + a). 


Thus, Uy — Do (x? + x?) > 0 if Do is sufficiently small, which gives (3.4) and hence 
by (3.2) and (1.8), leads to the estimate : 


V(t) > Dy (x? + x*) —{4g| x] |x| + Ar 2°} 


> Dy x* + (Dy — Ai|2) x2 — Ay |x| — 2A, G5) 
> Dy (x? + x) - $4, ¥*—D, 
for some D,, D, if A; is taken sufficiently small. 
Because of the (assumed) w-periodicity of x (t), we have, on integrating (3.5), that 


w w 


0>D, | (x? + x*)dt — 8 Ay| x? dt — D, w. (3.6) 
0 


0 
Combined with the inequality 


w 


| x? dt <4w? | x? dt es 
0 0 


which can be verified by substituting Fourier expansions of x and x in (3.6), (3.7) 
leads to the estimate 


w w 


(D, —§ w? 1-? Aj) [x dt + p,| x? dt < D, w. 


0 0 
Therefore, if A, is further fixed such that 
Aiw*n? < £D, 
as we assume henceforth, then 


w 


D, J (x? + x) dt < 2D, w. ...(3.8) 


In particular 


o- = 


0 dtie, (De 
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Considering now identity 


‘ 


x(t) = x(T;) + | x? (5) ds 


with 7, fixed (as is possible in view of the periodicity condition « (0) = *(w)) such 


that x (T,) = 0, we have that 


w w 


max | x(t)| < | | x(s)|ds wir (| x? (s)ds)'2 
0<t<w 0 0 


by Schwarz's inequality. Thus (3.8) implies that 


max | x(t)| < Di/? w. ...(3.9) 
0<t<w 


From this, on referring to the identity 


t 


x (t) = x (T») +| x (s) ds. 


With T, chosen such that *« (7,) = 0 (the choice being possible in view of the periodi- 
city condition x (0) = x (w)), we have that 


max | <(t) | & Di!? w?. .«-(3.10) 
0<t<w 
To obtain an estimate for | x(t) | first note that, because of the w-periodicity 
of x (t), integration of both sides of (3.1) yields the result 


w w 


J rey mnt rece x ai dime | fs (3) ai 


0 


or indeed, in view of (3.10), that 


| | ‘fie (x) ae Lp (t, x x, x, x)} dt <S D, ae a 


0 
by (1.8), (3.9) and (3.10) 
| Lp (t, x; x, x, x) | < Dit Ww 


for some D,. Also, since a, > 0 and f, (x) sgn x + 00 as | x | > 09 (by(1.7)), there 
clearly exists Ds independent of » such that 


| fasp (x) — PP (t, x, x, : x) | > 2D, wo -.-(3.12) 
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if | x(t) | > D, for all ¢ € [0, w]. It is thus clear that | x (T;) | < D, for some Ts, 
as otherwise, by (3.12), the left handside of (3.11) would be not less in magnitude than 
2D,. The result that | x (73) | < D; combined with (3.10) to yield the required 
boundedness estimate for x: 


max | x(t) < D, + Dj!’ w*. seis 35) 
O0<tSw 


It remains now to obtain estimate for | x(t) | in order to complete our verifica- 
tion of (2.2). Multiplying (3.1) by x“ and integrating from tf = 0 to tf = w we obtain 


w w w 


| oxenye die Gi | fi (x) xx de at x) Qdt ...(3.14) 
0 0 0 

where 
O = uplt, x x, x) — (fae (0) x — pfs (0) — fan (X)} 

satisfies 


|O|<D, (3.15) 


by (1.6), (1.8), (+.9), (3.10) and (3.13). But by (3.9), | fi (x) | < Dy, so that 


=n] AG) ex ar] <D, (f (oy any 


by (3.8). Thus from (3.14) and (3.15), we have that 
J (xl) dt < Da ( i (x9? ary? 

and hence that 
f (x5*)* dime" De 

for some Dy. The result 


max |x (t)| < w/? Dil? 
w 


aes ...(3,16) 


now follows readily. 


The result (3.9), (3.10), (3.13) and (3.16) full 
w-periodic solution x (t) of (3.1) if the A, in 
pletes the proof of the Theorem. 


y verify (2.2) for the arbitrary chosen 
(1.8) is sufficiently small, This now com- 
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If L (x, y) is the metric function of a Finsler space F" = (M*,L) and 8 is a 
one-form 5: (x) y* in F”, then the trannsformation L > L* = f (L, 8) known 
as &-change has been introduced by Shibata®. In the present paper we 
consider a particular case of 8 changes known as cubic transtormation given 
by L > L* = (L3 + 8°)!/3, The relation between n fundamental forms of 
tangent Riemannian hypersurface of (M",L) and (M",L*) has been 
obtained. 


1. INTRODUCTION 


Let F" = (M", L) be ann dimensional Finsler space with fundamental metric 
function Z (x,y). In general L (x, y) isa function of point x (= x‘) and element of 
support y (= y') and positively homogeneous of degree one in y. Let (M", L’) and 
(M", L") be Finsler spaces whose metric functions L’ (x, y) and L” (x, y) are obtained 
from Z by the relations 


Li'=L+8 : Sth 
bot ee lets ---(1.2) 


where B = 0; y‘, b; (x) isa component of a covariant vector which is a function of 
position alone. These two iransformations have been introduced by Matsumoto! 
which have the geometrical properties stated therein. Generalizing * these transforma- 
tions, Chibata® has studied the Properties of Finsler space (M*", L*) whose metric 
function L* (x, y) is obtained from L by the relation L* (x, y) = J (LZ, 8), where / is 
positively homogeneous of degree one in Z and 8. This change of metric function 


is called a B8-change. A particular case of 8-changes is a cubic transformation of 
metric function given by 


ee aes (1.3) 
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The m fundamental forms of a Riemannian hypersurface of a Riemannian space 
have been defined and their properties have been studied by Rund’’*®. Prasad® has 
obtained the relation in m fundamental forms of tangent Riemannian hypersurfaces of 
(M*, L) and (M’, L’) whose metric functions are related by (1.1). He® also obtained 
the relation in n fundamental forms of tangent Riemannian hypersurfaces of (M", L) 
and (M", L”) whose metric functions are related by (1.2). In this paper we shall 
obtain the relation inn fundamental forms of tangent Riemannian hypersurfaces of 
Fs = (M", L) and F*" = (M", L*) whose metric functions are related by (1.3). 


2. Tue RELATION BETWEEN v-CURVATURE TENSORS OF F" AND F*" 


Throughout this paper the quantities corresponding to the Finsler space F*" will 
be written by putting ‘*’. For the cubic transformation (1.3) the relation between 


the angular metric tensors /,;, and hy of F" and F**" will be given by 


he, = phi; + 2pq mm, aay) 


where p= LL*1,q =8 L*>,m = q I, — ph; and 1; = &L. Thus the metric 


tensor hi of F*" and its recipocal g*" are given by 


gt, = par + Pall, — pg? (hb; + 1 bi) 


+q(1 + p®) bb; net 2) 
gtd = p-gt — pgidA(p + gb) I'l + g(t! + Ib’) 
— 2pg » bt bd »20(2.3) 
where we put 
bt = gb, It = gl, b? = g% by by, 
h-) = p® — q® + 2p* gb’. 
The differentiation of (2.2) with respect to y* gives the (h) hy-torsion tensor 


C*, = (8 af, [a*)/2 of F*": 


ijk 


ce 


tik 


= pCijn + pq? (hij me + hie my + hy m)/2L 


— (2p? — 1) (p/L) mim, m. ..(2.4) 


Paying attention to m y‘ = 9, hi, y! = 0, Cijx y’ = 0, we obtain 


Cat Cg = ch, + (q?/2L) (hY m, + hy m) 


(equation continued on p. 244) 
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+ pq d (ql — 2pb') Cry + p? g® (p + 9b) (ZL) fu, PY 
— (pq? M2L) hy; b* — (qriL) [p® (2p* — I) (p + 95°) 


+ g®] mm, 1" + (pafL) mim, b" A Pe) 
where 
Ci; = C1 jk b*, 
We shall now find the y-curvature tensor rar which is given by S*,, = oe 
Crim — Cy Cin: By the use of (2.4) and (2.5) we get 
S¥ i, =P Sitje + Cry Ate + Cae diy — Cie day — Cry dix 
+ Ij Ene + Auk Eig — hix Eny — hry Eve ...(2.6) 
where 
dak = (p? AL) my me — pq? Cuz — (p? q? Al4LZ) hak Bal FE) 
En = (pq?aA/4L?*) [q? (2p? + 1) — 2p 5°] ma me 
— (pq! A[8L*) (g? — p* b*) hi — (p? q? A[4L) Crs. ira 


Remark : It is to be noted that if L (x, y) is the metric function of a Riemannian 
space i.e. L(x, vy) = (ai, (x) y‘ y’)'? = « then L* (x, y) = («® + B%)'/%. Thus F** 
is a Finsler space with an (a, 8) metric. It is known that a Finsler space with an 
«, B-metric is semi-C-reducible*. In order to clearify this fact in this case, calculating 
from (2.4) we have 


Coe pet (hy, Cr * hi, Cr = * hi Ci \(n + 1) 


VCE G Gels ...(2.9) 
where 
U = (n + 1) q? A/2L, V = — (1 + q*) C*? AL*, C*2 = gt CP CP, 


and 
A = 2L[{(n + 3 A) q? — 2p* Bb? A}. 


3. HyPERSURFACE OF (M", L) 
Let (M"~, L) be a hypersurface of (M", L) given by the equation 


xf = x‘ (u*). asCaabe 


Let us suppose that the functions (3.1) are at least of class C? in u* and the projection 
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factors Bi = dx‘/du* are such that their matrix has maximal rank n— 1. The 
fundamental metric function L (u, v) of the hypersurface is given by 


L (u, v*) = L (x* (ut), BL v*) 
where v™ is the element of support for the hypersurface for which 
yoo = Bi y*. ot 


Thus if /* denote the normalized vector along the element of support then 
Is = Bi I. If ga, (x, y) denote the metric tensor of (M"*, L), the induced metric 
tensor of (M"-1, L) is given by 


gap (u,v) = Bnj (X,Y) BL Bh. v..(3.3) 
The inverse of (3.3) is denoted by g*® (u, v) by means of which we define the quantities 

BS (u, v) = g®® (u, v) 81; (x, Y) Bs . ...(3.4) 

The unit normal vector N’ (u, v) of (/"~", L) is determined by the relations 

Bn, (x, ¥) BR N? (u, v) = 0, Bn, (x, ») N* (u, v) N?(u, ¥) = 1. 5) 
We have the following identities from (3.3), (3.4) and (3.5) : 

By; Bi = 65,5, B, + NIN = oH ...(3.6) 
where 


M: = 81; (x, y) N’. 


If Cu jx (x, ¥) denotes the (h) hy-torsion tensor of (M" L), the induced (h) hy- 
torsion tensor Caey (u, v) of (M*, L) is given by 


Capy (uy ¥) = Crye (X,Y) BL BE BY | ...3.7) 


from which we obtain 


e ge oad: ON: a 


. : t . 
The relative v-covariant derivative of the projection factor B, with respect to the 


induced Cartan connection /C TI is given by? 
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i « ] h k 
Bay =— B, Cp, + Cy 3, 4,- seat 


This tensor is normal to (M@"!, L). Therefore we may write 
aa i fat 
By, = May N’. (3.10) 
From (3.9) it is clear that Mgy is symmetric in 8 and y and it may be written as 


May = Ci N* By Be (ih ek 


The tangent vector space M?'' to M1 at every point x‘ (= u*) of the hyper- 
surface is considered as the Riemannian space (My" gx) with the Riemannian metric 
2x = Sak (u, vy) dv* dv®, The components of the (/:) hv-torsion tensor Ca, will be the 
Christoffel symbols associated with gx. If Mj isthe tangent vector space to M" at 
x‘ (= u*), then (M"", gx) will be the hypersurface of (M® , gx) given by (3.2), where 


&x = 81; (x, ») dy! dy’ is the Riemannian metric on M" . The quantities May given by 
(3.11) will be considered as the coefficients of the second fundamental form of the 


tangent Riemannian space (M"", 2x). 


_ In general the coefficients of the rth fundamental form of (M" , gx) are defined 
as" 


Cu )08 ="2aR; C208 = Mor, Cora = Cor 1)a8 Mg (2 = rs n) pil oet ae 
where 


Ma = gs Mup. 


4. THE n FUNDAMENTAL Forms OF A HYPERSURFACE oF (M", L*) 
Let (M"", L*) be a hypersurface of (M" L*) given by the same equation (3.1). 
It is to be noted that a unit normal vector N‘ to (M*""', L*) is not necessarily normal 
to(M"~?, L*). Paying attention to J, N‘ = 0, from (2.2) we have 
gt, By N* = (biN* {q(1 + p?) ba — p* q? Ia} (4.1) 


gi, N'NY = p+aq(l + p)(b, NY? (4.2) 
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where b, = 5; B.. If b, is tangential to the hypersurface (M"-}, 1), that is, 
b, N' = 0, then we have 

N*! = p-)? NE na 444) 
where we have chosen a normal vector N*‘ to (M""!, L*) in the same direction as N’. 


Hence we have the following : 


Theorem 4.1—Let (M", L*) be a Finsler space obtained from a Finsler space 
(M", L) by the transformation (1.3). If(M", L*) and (M"-!, L) are the hypersurface 
of these spaces given by the same equation (3.1) and 5; is tangential to the hypersurface 
(M*-}, L), then the vector normal to (M"-', L)is also normal to (M"~?, L) 


Now we establish the following : 


Theorem 4.2—Let (M", L*) be a Finsiler space obtained froma Finsler space 
(M", L*) by the transformation (1.3). Let(M"’,L*) and (M""',L) be the hyper- 
surfaces of (M", L*) and(M", L) given by the same equation (3.1). If db; is tangential 


to the hypersurface(M"-?, L), and (M", gx), (Mi g7),(M"”’, gy), (M”", gx), are the 


tangent Riemannian space to (M”, L), (M", L*) (M"-?, L) (M"™, L*) respectively, then 
we have the following : 


(i) The second fundamental forms of (M?* gx) and (M"~", g* )are proportional. 
(ii) Every asymptotic direction of Ce gx) is also an asymptotic direction of 


Hoole 9%): 


(iii) The rth fundamental tensors of (M""," gx) and (M"” g* ) are related by 


| cod | 
CoyaB = peo [Coyaa — Sy X (mB Qer+:—m)x] (3 Sr <2) (4.4) 
m=2 
where 
X(m)x = P /2ga Cimas 6 (2< man — 1) (4.5) 
Yom) = P /2gd Xims 0 (2 << m<n— 1) (4.56) 
Bead (Me 
Qryn = X(2)x (4.5c) 


r-} : 
Oryx = X(ryx —2% Ym) O(r4.1- mx (3 <rgan- I). (4.54) 
m=2 
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Proor: (i) Since b; N‘ = O implies m; N‘ = 0, from (2.4), (3.11) and (4.3) it 
follows that 


Mg, = p!* May. ...(4.6) 


This proves (i). 


(ii) A direction 1* for which Mag 1* 1° = 0 is said to be an asymptotic direction. 
In view of this definition and (4.6) we get (ii). 


(iii) The validity of relation (4.4) is established by induction. Since Cj, is an 
indicatory tensor, from (3.11) it follows that Mpy/” = 0. Hence from (3.12) we have 


CuryBy Iy=0 = Cu)py /8 (2:6 7s n). ...(4.7) 
Hence from (4.5) (a), (c) and (d) we get 
Xi rye 1* = 0, Oey. I* = O(2 |r < an). .-.(4.8) 


; e a pe 
Since g**® = g*!4 B® Bi, from (2.2) and (4.6) we get 

Mg’ = pl? [Mg + pq?d Map b* 18 — pr/2qgr Xp dF] ... (4.9) 
where b* = By b!. The relations (3.12), (4.5a) (4.6), (4.8) and (4.9) yield 


Cann = Conb— lg tas! ...(4.10) 


From (4.5c) and (4.10) it is evident that (4.4) holds for r = 3. Fora given fixed value 
of the integer s with 3 < s < n — 1 we have 


* & 
(s+1)0—B Clee Mg .. (4.11) 


Now let us suppose that (4.4) is valid for s = 3, 4, 5, ...r. so that we can write (4.11) 
in the form 


s~1 
x = 
Cit1)0B = p's)! (Cioas — SS X(m)8 Q541- mya] mM; 


m=2 
which in view of (4.5), (4.7), (4.8) and (4.9) gives 


scl 
bg = pl2- 5 
Cistt)aB sl “deh [Cis41)08 — X(m+1) Qs+1-m)o 
m=m2 


(equation continued on p. 249) 
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a=1 
— XB {Xia — 2 Yom Q(s41—m)a}] 


s 
= pl? [Coir08 — 2 Xeme Qs,2-mel.- 
m=2 


This shows that (4.4) is valid for r = (s + 1), which completes the proof of (iii). 


Remark | : If q¢ A < 0, the quantities X(,,)x, Yom) and Qy,,). defined in (4.5) are 
complex quantities. However the quantities Cr yop are always real. 


Remark 2: Theorem (4.1) and Theorem 4.2 are also valid for general 8-changes. 
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COPURE-INJECTIVE MODULES* 
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Department of Mathematics, University of Mysore, Mysore 570006, Karnataka 
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in this paper we introduce the notion of a copure-injective module and study 
its properties over special types of rings. Our main structure theorem asserts 
that over an arbitrary ring R, the copure-injective right R-modules are pre- 
cisely the direct summands of direct products of cofinitely related right R- 
modules. Over a commutative classical ring R the copure-injective R-modules 
are precisely the pure-injective R-modules. If R is a commutative ring then 
every pure-injective R-module is copure-injective. Over a Dedekind domain 
R, every copure homomorphic image of a copure-injective R-module is 
copure-injective. Finally we derive the analogue of Schanuel’s lemma for 
copure short exact sequences and copure-injective modules. 


Throughout this paper by a ring R we mean an associative ring with identity and 
by an R-module we mean an unitary right R-module while mod-R stands for the cate- 
gory of all right R-modules and R-homomorphisms. 


Definitions 1—(i) An R-module M is said to be ‘finitely embedded’ (Vamos"®, p. 
643) (or ‘cofinitely generated’ (Jans’, p. 588)) if E(M) = E(S,;) @E (S.) ® ... @ 
E(S,), where Si, S,, ..., S, are simple R-modules (here E (X) denotes the injective hull 
of an R-module X). 


(ii) An R-module M is said to be ‘cofree’ (Hiremath’, Definition 6) if M is iso- 
morphic to II {E (S.) : Sa is a simple R-module, « € A} for some index set A. 


(iii) An R-module 4 is said to be ‘cofinitely related’ (Hiremath’, Definition 14) 
if there is an exact sequence 0 > 4 > B > C > 0 of R-modules with B cofinitely 
generated, cofree and C cofinitely generated. 


(iv) A short exact sequence of R-modules is said to be ‘copure’ (Hiremath®, 
Definition 3) if every cofinitely related R-module is injective relative to this sequence. 


(v) A ring R is said to be ‘right co-noetherian’ (Jans, p. 588) if every homo- 
morphic image of a cofinitely generated R-module is cofinitely generated. 


(vi) A short exact sequence 0 > 4 > B>C—>0 of R-modules is said to be 
‘pure’ (Warfield®!, p. 703) if for every left R-module M, the induced map f @ Iu: 
A @rM — B @rM of abelian groups is a monomorphism. 


*This work is based on the author’s doctoral dissertation and was supported by the University 
Grants Commission and Madurai Kamaraj University. 
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| (vii) An R-module is said to be ‘pure-injective’ (Warfield”!, p. 703) if it is injec- 
tive relative to each pure short exact sequence of R-modules. 


Definition 2—An R-module is said to be copure-injective if it is injective relative 
to every copure short exact sequence of R-modules. 


From the definition we have following easy consequences, 
Proposition 3—(i) Every injective R-module is copure-injective. 
(ii) Every cofinitely related R-module is copure-injective. 


(iii) A copure short exact sequence 0> 4 > B>C->0 of R- 
modules with A copure-injective, splits. 


(iv) Every direct summand of a copure-injective R-module is co- 
pure-injective. 


(v) Every direct product of copure-injective R-modules is copure- 
injective. 


Remark 4: A copure-injective R-module need not be injective. 


Example—Let G = Z (n) be a cyclic group of order n. Since G is finite it is co- 
finitely generated and hence, cofinitely related as a Z-module (Hiremath’, p, 5 and Pro- 
position 17). So, by Proposition 3 (ii), G is copure-injective as a Z-module. But we 
know that G is not injective as a Z-module. 


Proposition 5—For a ring R the following conditions are equivalent: 
(i) Risa right V-ring; 

(ii) every cofinitely related R-module is injective; 

(iii) every short exact sequence of R-modules is copure; 

(iv) every copure-injective R-modules is injective. 


Pooor : (i) <> (ii) < (iii) follow by Proposition 5 of Hiremath®. (iii) > (iv) is 
obvious and (iv) => (ii) follows by Proposition 3 (ii). 

Remark : A direct suin of (copure—) injective R-modules need not be copure- 
injective (see example (i) below). We know that over a right noetherian ring R every 
direct sum of injective R-modules is injective (Faith*, Proposition 6.5). But over a 
right noetherian ring R (even when & is a Dedekind domain) a direct sum of copure- 
injective R-modules need not be copure-injective (see example (ii)). 


Examples—(i) Let R = MM Z/p, Z, where {p;, Py, .... Pn» +} is the set of al 
a=] 
positive primes. Since Ris a Von Neumann regular ring, R isa V-ring by Theorem 


6 of Rosenberg and Zelinsky'®. Let, form = 1,2, .., @, = te Zip, Z. Then Q; a 
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Qc... C Gs C «» isa strictly ascending chain of ideals of R. Then, by Proposi- 
+ + + 


io 2) 
tion 6.5 of Faith’, E = ® E(R/Q,) is not injective and so, not copure-injective by 
n=1 


Proposition 5, as R is a V-ring. 


oo 
(ii) LetG = ® Z/p, Z where {Pi, Ps, ..-, Pn, --.} is the set of all positive 
n=l 


primes. Since each Z/p, Z is finite, it is cofinitely related as a Z-module by the remark 
on p. 5 and Proposition 17 of Hiremath’ and so, by Proposition 3 (ii), each Z|PyrZ 18 
copure-injective as a Z-module. Clearly G isa pure submodule of the Z-module 


lo e) 
H = 1 Z/p,Z. Hence, by Proposition 12 of Hiremath®, G is a copure submodule 
n=l 


of the Z-module H as Z is a neotherian ring. If G were copure-injective, then by Pro- 
position 3 (iii), G would be a direct summand of H whichis absurd. SoG is not 


copure-injective as a Z-module. 


We recall (Faith*, p. 254) that a category @ is ‘locally small’ if the equivalence 
class of subobjects of any object of @ is a set and @ is ‘colocally small’ if the dual 
category @* of @ is locally small. Since the category mod-R is colocally small (Faith*, 
Exercise 5.27.4), the isomorphisms closed class of cyclic R-modules, distinct up to iso- 
morphisms, is a set. Hence the isomorphism closed class ¢ of simple R-modules, 
distinct up to isomorphisms, is a set. Since mod-R is a locally small category (Faith‘, 
Exercise 5.27.4), the isomorphism closed class of subobjects of the R-module A = 
Tm {E(S): S € &} (distinct upto isomorphisms) is a set. Since any cofinitely generated 
(and hence, cofinitely related) R-module is isomorphic to submodule of A, the isomor- 


phim closed class C* of cofinitely related R-modules, distinct up to isomorphisms, is 
a set. 


We now prove that for any ring 2, there are enough copure-injective R-modules 
in the sense that every R-module can be embedded in a copure-injective R-module as 
a copure submodule by applying Proposition 3 to: ; 


Theorem 7—Every R-module can be embedded ais a copure submodule in a 
direct product of cofinitely related R-modules. 


Proor : Let C2 be the isomorphism closed class of all cofinitely related R- 


modules distinct up to isomorphisms. We have seen above that C* is a set 
r . 


Let A be any R-module and let 


Hom p(4»M) 


Q=I1{M : ME C*}. 
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Define a map 


¢:A—>Qby 
$ (a) = ((f (4) feHome(4.My) ME C” a EA. 


Clearly ¢ is an R-homomorphism. We claim that ¢ is a copure monomorphism. Let 
OAa EA. Let D7be a maximal right ideal of R containing annr(a). Then 


S = R/D is cofinitely generated and hence E(S) € C* . Define a map«: aR > S 


by «(ar)=r+ D7,reER. Then«x extends to a homomorphism f: A — EF (S). 
Since f (a2) = 1+ DT £0, ¢ (a) #0, proving that ¢ is a monomorphism. 


To prove that ¢ is copure, let M € Os and f © Homer (A, M). Then clearly 


f = pd where p: O>Misthecomposition zy o my of the usual projections 7m and my. 
We now have the following structure theorem for copure-injective modules. 


Theorem 8—An R-module M is copure-injective if and only if it is a direct sum- 
mand of a direct product of cofinitely related R-modules. 


Proor : Necessity follows from Theorem 7 and Proposition 3 (iii) and sufficiency 
from (ii), (iv) and (v) of Proposition 3. 


We now compare copure-injectivity with pure-injectivity. For this we first 
observe that from the adjoint isomorphism of Hom and @ one can easily deduce that 
for a left R-module M, the canonical R-module M* = Homz (M, Q/Z) is pure-injec- 
tive. 


Proposition 9—Over a commutative ring R every pure-injective R-module is co- 
pure-injective. 


Proor: This follows from the fact that fora commutative ring R, copurity 
implies purity (Hiremath®, Proposition 13). 


Remark 10: A pure--injective R-module need not be copure-injective. 


Example—Cozzens? has constructed a ring R = k [x, D] of differentiable polyno- 
mials in a single indeterminate x over a universal field k with a derivation D (where 
multiplication is given by ax = xa + D(a), a € k). Cozzens has proved that Ris a 
right V-ring but not a field. Since x is not invertible in R, Rx is not a pure left ideal 
of R and hence R/Rx is not flat as a left R-module. So by the Corollary on p. 131 of 
Lambek!”, (R/Rx)* = Homz (R/Rx, Q/Z) is not injective as an R-module whence, by 
Proposition 5, (R/Rx)* is not copure-injective. But (R/Rx)* is pure-injective. 


Proposition 11—If R is a commutative (co-) noetherian ring then every copure- 


injective R-module is pure-injective. 
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Proor : Since over a commutative (co-) noetherian ring R, every pure submodule 
of an R-module is copure (Hiremath®, Proposition 12), the proposition follows. 


We now prove a proposition which dualizes (i) <> (iv) of Proposition 5: 


Proposition 12—A ring R is a Von Neumann regular if and only if every pure- 
injective R-module is injective. 

Proor : Since for a Von Neumann regular ring R, every short exact sequence of 
R-modules is pure, by Theorem 11.24 of Faith‘ the ‘only if part follows. 


For the ‘if’ part, suppose that every pure-injective R-module is injective. To prove 
that R is Von Neumann regular, we need only prove, by Theorem 11.24 of Faith‘, that 
every left R-module M is flat. Indeed, since the R-module M* = Homz (M, Q/Z) is 
pure-injective, it is injective by hypothesis. It then follows. from the Corollory on 
on p. 131 of Lambek'*® that Mis flat asaleft R-module proving that Ris Von 
Neumann regular. 


Remark : A copure-injective R«module need not be pure-injective. 


Example—Let V be a countably infinite dimensional vector space over the field 
Q of rational numbers. Let R be the ting of linear operators of V. Then R isa Von 
Neumann regular ring. Let PZ be a maximal right ideal of R containing the two-sided 
ideal J of R of all elements of R ot finite rank, Then, by Theorem 1 of Osofsky, 
S= R/DH is not injective as an R-module. Let xEE (S)\\S and let A be a submodule 
of E(S) maximal with respect to S C A and x @& A. Then E(S)/A is subdirectly irre- 
ducible and hence, is cofinitely generated. So Ais cofinitely related and by Proposition 
3 (ii), A is copure-injective. Since A is not a direct summand of E (S), A is not injec- 
tive and so not pure-injective by Proposition 12. 


From Propositions 9 and 11 we have: 


Corollary 14—If R is commutative (co-) noetherian ring then R-module is pure- 
injective if and only if it is copure-injective. 


We now recall the following Definitions. 


> 


(i) (Warfied*’, p. 707-708). Let J be any index set and let M be any R-module. 
Let M! be the /th cartesian power of M. For any finite subset J* of J, and elements 
ri€ Ri € 1*), we define a group homomorphism (which will be a module homo- 
morphism if R is commutive) 

@: M' - M by ¢ (x) = = x, Bya linear equation over Rin M we mean 


1 
a pair (4, m) where ¢ is a homomorphism of the type defined above andm € M. The 
solution set of this equation S (%, m) is the set of all elements x € M' such that ¢ (x) 
= mM, 
A family of linear equations over R in M is said to be finitely soluble if the cor- 
res p onding sets S (¢, m) have finite intersection Property. An R-module M is said to 
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be algebraically compact if every finitely soluble system of linear equations over R in 
M has a simultaneous solution. 


(ii) An R-module M is said to be linearly compact (in the discrete topology) (Warfield?!, 
p. 711 and Vamos”, p. 115) if every family of cosets in M with finite intersection pro- 
perty has nonempty intersection. 

(iii) A commutative ring R is said to be classical (Vamos*®, p. 121) if E (S) is 
linearly compact forevery simple R-module S (or equivalently, every cofinitely 
generated R-module is linearly compact). 


We now generalize the Proposition 11 to commutative classical rings by noting 
that commutative (co-) noetherian rings are classical (Vamos'*, Theorem 2 and Vamos?°, 
Proposition 4.1). 

Proposition 15— Over a commutative classical ring R, every copure injective 
R-module is pure-injective. 


Proor : This follows from Proposition 9 of Warfield”) that over a commutative 
ring R, every linearly compact R-module is algebraically compact and hence, pure-in- 
jective (Warfield*’, Theorem 2), from Theorem 8 and the facts that pure-injective R- 
modules are closed under taking arbitrary direct products and direct summands. 


Corollary 16—For a commutative classical ring R, every pure shart exact sequ- 
ence of R-modules is copure. 

Corollary 17—For a commutative classical ring R. an R-module is pure-injective 
if and only if it is copure-injective. 

Proor: This follows from Propositions 9 and 15. 


We next characterize the rings R for which every R-module is copure- 
injective. 
From Theorem 8 we have : 
Corollary 18—For a ring R the following conditions are equivalent. 
(i) Every R-module is copure-injective. 
(ii) Every R-module is a direct summand of a direct product of cofinitely related R- 
modules. 
(iii) Every copure short exact sequence of R-modules splits. 


Fieldhouse® (p. 15) calls a ring Ra ‘right PDS ring’ if every pure submodule 
of an R-module is a direct summand (or, equivalently (Fieldhouse®, Theorem 10.1) 
every pure short exact sequence of R-modules splits). Similarly we call a ring a right 
CDS ring if R satisfies the equivalent conditions of Corollary 18. 


Since a ring R is a right PDS ring if and only if every R-module is pure-injective, 
we have from Proposition 9 and Corollary !7: 
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Proposition 19—(i) If R is a commutative PDS ring then R is a CDS ring, 


(ii) If Ris a commutative classical ring then Risa PDS ring 
if and only if it is a CDS ring. 


We now have, from Proposition 19 and Theorem 10,3 of Fieldhouse : 
Corollary 20—Every commutative classical CDS ring is artinian. 


We do not know whether a commutative artinian ring is CDS. However, we have 
from Theorem 10.4 of Fieldhouse : 


Corollary 21—Every commutative uniserial ring, that is, a commutative artinian 
principal ideal ring, is a CDS ring. 


Remark 22: We know that Dedekind domain need not be a CDS (= PDS) ring 
(cf. Proposition 19 (ii)) (e. g. Z, the ring of integers, is nota CDS ring since Z is 
not algebraically compact as an abelian group (Fuchs®, Chapter VII, § 38, Exercise 1)). 
But every proper homomorphic image of a Dedekind domain, being an artinian prin- 
cipal ideal ring (and hence a uniserial ring), is a CDS ring by Corollary 21. 


Remark 23: We know (Cartan and Eilenberg’, Proposition 6.!) that a ring R 
is right hereditary if and only if every homonorphic image of an injective R-module is 
injective. But neither a homomorphic image of a copure-injective R-module over a 
right hereditary ring need be copure-injective (see example (i) below) nor a ring R for 
which every every homomorphic image of a copure-injective R-module is copure-injec- 
tive need be right hereditary (see example (ii)). 


Examples—(i) We know that Zisa hereditary ring and copure-injective Z- 
modules are precisely the algebraically compact Z-modules (Theorem 2 of Warfield?! 
and Corollary 14). Since there exists an algebraically compact abelian group with a 
homomorphic image which is not algebraically compact, it follows that there exists a 
copure-injective Z-module with a homomorphic image which is not copure-injective. 
(ii) R = Z/(4) is an artinian principal ideal ring and hence is a uniserial ring. Then 
it is a CDS ring by Corollary 21. So, every R-module, in particular every homomor- 
phic image of a copure-injective R-module, is copure-injective. But R is not 
a hereditary ring since the maximal ideal (2)/(4) of the local ring R is not projective as 
it is not free as an R-module. 


Remark 24: We know (Fuchs®, §38, Exercise 3) that every pure homomorphic 
image of an algebraically compact abelian group is algebraically compact. But, in 


general, a copure homomorphic image of a copure-injective R-module need not be 
copure-injective, 


Example—Let R = I Z/pZ where P is the set of all positive primes. Since R 


. . . . . P . . 
is self-injective (Faith, Exercise 5.64.2), it is self-copure-injective. Since R is not 
semisimple artinian, there exists, by the Theorem of Osofsky"’, an ideal J of R such 
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that R// is not injective as an R-module. Since Risa V-ring, R/I is not copure-injec- 
tive as an R-module by Proposition 5. 


We prove below, in Corollary 28, that every copure-homomorphic image of a 
copure-injective R-module over a Dedekind domain R is again copure-injective. Before 
proving this we take note of and derive a few results: 


(i) Every finitely presented R-module is pure- projective. 


(ii) Any R-module can be embedded in a pure short exact sequence 0 — K -> P > 
M — 0 of R-modules where P is pure-projective. 


(iii) An R-module is pure-projective if and only if it is a direct summand of a direct 
sum of finitely presented R-modules. 


Next we recall (Kaplansky”, p. 332) that if R is a commutative integral domain, 
then an R-module M is said to be decomposable if Misa direct sum of cyclic R- 
modules and finitely generated R-modules of rank one. Since over a Dedekind domain 
R, every finitely generated torsion-free R-module is projective (Cartan and Eilenberg’, 
Chapter VII, Proposition 4.1 and a remark on p. 134), it follows that every decom- 
posable R-module, over a Dedekind domain R, is a direct sum of cyclic R-modules 
and a projective R-module and hence it is pure-projective by (i) and (ii) of the results 
stated above. 


Proposition 25—Over a Dedekind domain R every submodule of pure-projective 
R-module is pure-projective. 


Proor : We first observe that every pure-projective R-module is decomposable. 
This follows from: (a) Steinitz’s theorem (Steinitz'”) that every finitely generated 
R-module over a Dedekind domain R is decomposable; (b) the decomposability is a 
hereditary property for a Dedekind domain (Kaplansky’, Theorem 4) and (c) form 
(iii) of the properties of pure-projective modules stated above. 


From this observation and the observation made just before this proposition it 
follows that pure-projectivity and decomposability are equivalent for a Dedekind do- 
main. The proposition now follows from Theorem 4 of Kaplansky’?. 

Proposition 26—An R-module Q is pure-injective if and only if O is injective 
relative to each pure short exact sequence 0 > A>B>+C+0 of R-modules with 
B pure-projective. 

Proor: This follows by using arguments similar to the Proposition 5.2 of Chapter 
] of Cartan and Eilenberg! and the fact that the class of all pure short exact sequences 
of R-modules form a ‘Proper Class’ in the sence of MacLane’® (p. 367) (Stenstrém'’, 


Propositions 2.2 and 9.1). 
By Propositions 25 and 76 and using arguments similar to Theorem 5.4 of 
Chapter J of Cartan and Eilenberg!, we have: 
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Proposition 27—Every pure homomorphic image of a pure-injective R-module 
over a Dedekind domain R is pure-injective. 


Since, by Theorem 20 of Hiremath®, purity and copurity are equivalent fora 
Dedekind domain we have, from Corollary 14 and Proposition 27 : 

Corollary 28—For a Dedekind domain R, every copure homomorphic image of a 
copure-injective R-module is copure-injective. 


Finally we prove the analogue of Schanuel’s lemma for copure short exact seque- 
nces and copure-injective modules. We give the proof for the sake of completeness. 


i g 
Proposition 22—Let A be an R-module and let 0 > A +> O +> B>Oand0>4A 


Ed Q' ce B -> 0 be copure short exact sequences of R-modules with Q, QO’ copure- 
injective. Then Q @ B’ and Q’ @ Bare isomorphic. 


ProorF : Considering the usual pushout C =(Q@ O’)|\K where K = 
{(f (@),—f' (a)) : aE A} and the natural maps fh, h’ and by defining k, k’ in the obvious 
way we obtain the following commutative diagram of exact rows and columns : 


0 
| fl g 
0— A— O— B—> 0 
pla] « |b 
lm 5 Fn Co ote eg) 
“| Ip | F 
B’'— > B' 
0 0 


We prove that the Second row and the second column are copure. Let M be any 
cofinitely related R-module and let « : Q’ — M be any homomorphism. Then by the 
copurity of the first row there is a B: QO + M such that Bf = af’. Defined:C > My 
by ¢ ((x, x‘) ++ K) = 8 (x) +4 (x’) for (x, x’) + Kin C. Clearly ¢isa well-defined 
homomorphism and ¢ h’ = «. So the second row is copure. Similarly the copurity of 
of the second column follows. 


Now by the copure-injectivity of O, O’ and the copurity of the second row and the 
second column, both the second row and the second column split by Proposition 3 (iil). 
Hence C is isomorphic to both Q ® BandQ’ OB. 
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SEQUENCE SPACES OVER NON-ARCHIMEDIAN FIELDS IV 
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Necessary and sufficient conditions for an infinite matrix defined over a 
field K with non-trivial non-archimedian valuation to transform I* (K) into 
Co (K) andT (K) to I'* (K) are investigated. Tne sequence spaces Co (K), 
T (K) and I’* (K) are all defined over such a field K. 


1. INTRODUCTION 


The object of the present paper is to obtain some inclusion theorems of certain 
sequence spaces over non-archimedian fields which were not considered earlier?-®. In 
section 2 we shall describe the required preliminaries, whereas in section 4, we shall 
prove the main theorem of this paper. 


2. ENTIRE FUNCTIONS Over K 


Let K be a non-archimedian non-trivially valued field which is complete under 
the metric of real valuation. Let Ny be defined as 


N={|x|:x€ RK}. 


Leta (x)=a+a,x+ax+...+a,x"+... bea power series with coefficients 
in K. As noted by Raghunathan? (p. 517), « (x) is an entire function over K if and 
Only if 

| a, |?/" + 0 asn > 0, 
If « (x) is an entire function, let 


M (R) > aie | a (x)|,(R © M) 


Then by a known result! (p. 85), we have 
. M (R&) = “n= 
( so | a, | R",n > 0] ---(2.1) 


where 


ic,2) 
oo >} a, x*. 
n=0 
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3. DEFINITION OF DIFFERENT SEQUENCER SPACES 


In what follows, the notion of convergence and boundedness will be in relation 
to the metric of valuation of the field. 


C,(K) : The set ofall null sequences x = (x,) 
I'(K) : The set of all sequences x = (x,) with 
|x, | 1/" > O0asn > o9. 
['*(K) : The set of all sequences x = (x,) with 
{ | x, | ?/"} bounded. 
If x = (x,) is a sequence over K, let us define || x || = a | x, |. The norm is 
n 


evidently non-archimedian in the sense that it satisfies stronger form of triangle in- 
equality. With this as the norm C (K) is a non-archimedian Banach space. 


Let x = (x,) © K. Let | x | = Sup{|x,|1,.2 > 1}. Then | x | satisfies the 
following conditions. 


(i) |x| > 0, |x| = Oif and only if x = (0,0, ...... ) where 0 is the zero 
element of the field. 


(ii) |x +y |< Max{|x|, |yl} 
(iii) |tx| < A(t)|x|,¢€ K, A(t) = Max {1, | ¢| 3. 


Hence | x —y | defines a metric in the set of all sequences I (K). Sol (K)is a metric 
space defined over K, the topology being defined by the metric given above. 
Raghunathan? (p. 518) has proved that I (K) isa complete linear metric space which 
is totally disconnected. 


According to Raghunathan? every continuous linear functional f (x) defined for 
x€ET (K) is of the form f(x) = = C,Xn, X = (x,) where { | C, | *!"} is bounded. Hence 
I'* (K) is identified as the dual space of '(K). It was known from Raghunathan? 
(p. 524), I'* (K) is a complete metric space which is not linear. 


4. MATRIX TRANSFORMATION OF I'* (K) INTO Co (K) 
Let us consider the matrix transformation 
y= 3 ao Xp, 0 = 1, 2, 3, ..., and Opp E K, ..(4.1) 
pol 


Theorem 1\— A necessary and sufficient condition that (y,) © C, (K) whenever 
(x,) € T* (K) is that 
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{f, (x)} is bounded uniformly on every finite circle 


|xl<R,RE Me ...(4.2) 


where 
co 
tn (X) — >» anp Xpy A= Ie fs 3 coe and Anp E K 
p=l1 


is a sequence of entire functions over K. 
Proor : The condition is sufficient. 


Let 


fr (x) = = np Xp, Nn = 1, 2,3, ... anda,p € K 
p-~l 


be uniformly bounded on every finite circle | x | < R, R € Nz. Then the condition 
(4.2) implies that there exists a constant M (R) such that | f, (x) | < M (R) for every 


x with the property that | x | < R,R € M,. By using (2.1) and the condition (4.2), 
we have 


| Gnp | R? < M (R) for each fixed p. ...(4.3) 
Since (x,) € I'* (K), there is a constant c such that 

| x» 17!" < eso that | x, | < c? forall p. ...(4.4) 
From (4.1), 


[yn] Sup [ay] | xp}. 
1S poo 


By using (4.3) and (4.4) in the above, we have 


M(R) 
Re : 


Choose R = ye where & E Neand p > 1, 


I\ml< 


- 


M (uc) M (uc) 
Then | y, | < (H Cc? = =the) < « for large p so that y, > Oasn > co, 


The condition is necessary. 


Since fn (x) = % dnp x? is defined for every n and x = R, the series is convergent 
for [x= R: 


| Gap 2? | = | a,p | R? +O as p > o0, 


-..(4.5) 


When we consider the special sequence (0, 0, ..., 0, 1, 0, ...) where 1 is in the 


pth place, 


Yn = Gap, n = |, B 3, wren 
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Since (y,) € Cy (K), we have for each fixed p 


an,p > Oasn > co, 


.. (4.6) 


If the condition (4.2) is not satisfied, there is some finite circle | x | <R, REN; 
such that { f, (x)} is not bounded uniformly on | x | = R. Since | x | = Ris bounded 
in the metric of valuation, we can find a sequence (x,) such that for | x, | = R 
Tn (%n) > e° as n + o9, through a subsequence of values of 7. 


Using this assumption, we shall construct a sequence (x,) € I'* (R) with the 
additional condition | x, | = R’ for every p such that the corresponding (y,) is not 
bounded. This proves the necessity of the condition. 


Choosing n, such that 
| ny (Xn) | > € for some e > 0. .. (4.7) 


Since | x,, | = R, (2.1) and (4.7) together imply, 


Sup | a, ,|R? > «. .. (4.8) 
Ixp<o 


By using (4.5), we have 


Sup | dp | RP < 3 for D > Pn: (4.9) 
Pn, t1<p<e 


From (4.8) and (4.9), we get 


Sup | nip | RP > €. 
1SPSPn, 


Hence there is a p, in 1 < p < P,, such that 
1G aap Pd doom Oe ...(4.10) 


Now choose x, in] <p < P,, as follows. 


if forp=p,and[x|=R AAD 
xp = ; 
f 0 for all pin! < p< pm. 
Then 
a 
co 
Le Bow Z AnpXp + YB An, Xp. epee 


€ 
| i Qn ,p Xp l< Sup | Qn,» Il Xp Wo > by (4.9). ...(4.13) 
Pasa Pn tisp<% 
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Hence we shall have from (4.12) 


PoP. 


1 co 
| 2 A, pXp| = | np Xp |< Max{ly, 1,1 % An p Xp | }- 
p=l eo el 1 1 ao 1 


By (4.10), (4.11) and (4.13), we have 


« < Max {1 ym |») 
Therefore from the above, we have 
LePead es 
Since 4, — 0asn— oo for each fixed p. So | a,, | R? > 0asn — co 


for each fixed p 
| a,p | R? < €foralln > ny. 


Now choose 1, > m such that 
Sup |4,,p | R? > € } 


IS p< 
and 


| 
| 
i r 
Sup? PG. 1 a 
ISPSPn, 2 J 
This is possible if n, is large enough such that 
nm, > Max (n,) where 1 < pC p,,. 


Then there exists by (4.5) p,. greater than Pn, Such that 


€ 
Sup | Quop | R?< =, 
Pn, t1<p< 2 


Therefore from (4.14) and (4.15), we have ; 


SUp*, fa, (Re Se 
ISpSpny : 


Hence there exists a ps in 1 < p < p,. such that 
| nop, | Ree 4, 


So p, chosen in (4,17) exceeds Pn, by (4.14). Now choose x, as follows. 


7 Fae = R 
ed 
(0 for all pin p,, + 1 SPS Puy 


..(4.14) 


. (4.15) 


...(4.16) 


-+-(4.17) 


(4.18) 
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p“p 
ny =f 
_ 2) 
Yn, = & Qn p Xp —_ p> Anop Xp + x An p XP 
pal 2 +1 p +1 
ny no 2 
Pp 
no P, 
2 ag pXp = z 
nP ei np tp — xz a Xp. 
or 2 pm 2 > +) Hgts 
n2 2 
By using (4.18) in the above, we get 
my, 
>> = 
| . Paint Xp | | Any Po Xpe | 
RI 
Pp 


nl 
: co 
<= Max {e] Ye, |. | 2 Gp Xe |» | 2 Anop Xp | } 
t=) = fe +1 


By (4.15), we have .. (4.19) 
= < 
1g Se hi a aaa (4.20) 
i 2 
By (4.14), we have 
Pp 
™ 
|= a,,xXp| < Sup la, >| R°< —. 
aes I<p<pr 2 
1 
Using (4.17), (4.18), (4.20), (4.21) in (4.19), we get 
€ < Max erase at 
This implies | y,, | >€. 
e. That is 


Proceeding in the same manner, we can find n, such that | yn, | > 


| Yn, | does not tend to zero as n > oo through a subsequence of values of n. This 


shows that (y,) does not belong to C, (K), while the sequence (xp)E I'* (K). This con- 
tradiction proves the necessity of the condition. 


By using a method similar to the proof of the above Theorem 1, we can prove 


following theorem. 
Theorem 2— A necessary and sufficient condition that {y,} should belong to 
l (K) whenever (x,) belong to I'* (K) is that | f, (x) | 7" > Oasn > 0° uniformly on 


266 


D. SOMASUNDARAM 


every finite circle | x | < R, R © Nx where 


tin (x) = S Gap x? x nen 293, ae 
pai 


is a sequence of entire functions over K. 


Se Niet Man 
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COMMENCEMENT OF COUETTE FLOW IN OLDROYD LIQUID WITH 
HEAT SOURCES 


G. C. Dasu 


Department of Mathematics, College of Basic Science and Humanities 
O.U. A.T., Bhubaneswar 751003 


AND 


S. BiswAL . 
Department of Physics, Nayagarh College, Nayagrh 752069, Puri 


(Received 11 December 1987; after revision 14 June 1988) 


Commencement of Couette flowin Oldroyd liquid with heat sources has 
been analysed. Approximate solutions are derived by Galerkin technique. 
The effects of various parameters on the velocity field, the temperature field, 
the skin friction and the rate of heat transfer are discussed. It is observed 
that the flow is sensitive to the interactions of heat source and elasticity of 
the fluid. 


1. INTRODUCTION 


The problem of Couette flow between parallel walls has a rich stock of literature. 
Pai’ has discussed the problem of unsteady Couette flow of a viscous inrompressible 
liquid between two plate walls occurring due to the sudden motion of one of the walls. 
The technique used for the solution of the equations is the principle of superposition 
and the solution occurs in the form of an infinite series of functions of the variables 
Y and t. Nanda* has studied the same problem assuming the walls to be porous. 
Then, it has been extended by Katagiri? and Muhuri* Rath et al.5 who have solved 
the heat transfer problem in case of unsteady Couette flow between two parallel walls 
having different temperature. 


Unsteady Couette flowin non-Newtonian fluids has enormous applications in 
technology. Therefore, many researchers have obtained the approximate solutions 
for the unsteady Couette flow caused in visco-elastic fluids with different wall velocities 
under different physical situations, Generalized plane Couette flow of one Oldroyd 
fluid with either suction or injection at the stationary wall has been discussed by 
Mishra®. Soundalgekar’ has analysed the problem of plane Couette flow of Walters 
liquid B’ with equal rate of injection at one wall and suction at the other (moving 
wall). Padhy® has analysed the commencement of unsteady Couette flow in second 
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order liquid. The object of the present paper is to study the commencement of 
unsteady Couette flow in case of Oldroyd liquid contained between two porous walls 
in the presence of heat sources. 


2. Basic EQUATIONS 
The visco-elastic fluid model considered here is the simplified form of Oldroyd’s 
B-liquid® whose stress-strain relation is given by 


) 
[1 + Ay | pitk = 2% [1 + A, Sr ef; tL) 


In this equation A, and A, are the time of relaxation and retardation respectively, mo is 
the co-efficient of viscosity, P’ and e“* are respectively the stress and the strain rate 
tensors, 5/d¢ denotes the convective derivative and, for a contravariant second order 
tensor, b‘* is given by 


dbtk ab‘* 





uu os m ik 
ot or “a Dom 


+ ee bre + vi btm, 
v‘ being the velocity vector and 

21k = Visk + Viste 
Neglecting the second order terms in A, and A: in eqn. (1) and following Gieskus’’, 
we get 

P't® a= Qn, e* — aes = et* ae) 
where 

Ky = No (Ar —. Ag). 
The momentum equations are given by 


[ vt 
p | ar se a vs, | mos ee = ps --- (3) 


- 


where fp is the density of the medium and pan arbitrary isotropic pressure. The 
equation of continuity is 


U;, =0 (4) 
and the energy equation is given by 
DT : 
pC, Di eee T+ ectoe 


where C, is the specific heat at constant pressure, k* the thermal conductivity and ® 
the dissipation function prescribed by 


ae ...(6) 
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3. FORMULATION OF THE PROBLEM 


Here, the unsteady Couette flow of an incompressible viscoelastic fluid begins 
with the sudden motion of the lower wall with time varying velocity At’ where n is 
positive. The x’-axis is chosen along the lower wall and the y’-axis normal to it. The 
upper plane is specified by the equation y’ = L where the number Z will be defined 
later. It is also supposed that the walls extend to infinity in both sides of the x’-axis. 
We consider the suction and injection velocity V at the walls to be constant. Then 
the velocity components u’ and v’ at any point (x’, y) in the flow field compatible with 
the equation of continuity are given by 

u’ = u' (y, t) } 


aD. (7) 


Under these conditions and also due to the shearing action of the fluid layers, the flow 
field and the temperature field with viscosity dissipation for an Oldroyd fluid are 
characterized by the following equations : 


ou’ ou’ 07u' 031,’ 034’ 
ep (a + VS )= tae — he (sora + V ar | cate) 
06’ BON Pel gt Or: ou’ \? e*u’ ou’ 
PCs (ae + hay) = mae ae) Serie oe 
du’ 07" ee 
oe aa : =r | + S’ (6’ — 42). «+(9) 


The notations P, Cp, yo, Ky, 6’ and K* used in the above equations are respectively the 
density, specific heat, co-efficient of viscosity, coefficient of elasticity, temperature at 
any point and thermal conductivity of the fluid. The last term in eqn. (9) represents 
the source-sink term with S’ as its strength. 


The adequate boundary conditions are 
t' = 0:u' = 0 forall y’ : 
u’ = At™ for y’ = a 


fos: 
u=0 fory=L ...(10) 
and 
t' =0: 6 =0 forall y’ } 
e6' 
—_ = ‘= l 
ay’ 0 for y 0 (1) 
ite a) Se 
6’ = 6, fory’ = L. J 
The fact ore = 0 when y’ = 0 implies that the lower wall is a non conducting 
6 


y . . . ~ 
one. Now, we introduce the following dimensionless variables and parameters : 
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y’ ra u’ } 

Aube im ne eer | 

R= oa Re = a ti Sanat ee lle) 
open pert 4S’ 

E Sdeai > PO ere | 

= ’ and L = vv, T | 


where TJ is some reference time, #6, the temperature of the upper plate, and R, R, E, 5 
and S are respectively the suction parameter, elastic parameter, Eckert number, Prandtl 
number non-dimensional source-parameter. 


By the use of (12), eqns. (8) and (9) are reduced to non-dimensional form as 
follows : 





eu e*u Ou c8y Fu 
e6 Cé@ 1 3%g O'u Cu Cu | 
eh gery ee Ha t ReEo[R patarid ast see 
Cu \? l 
A a Oe a 
+ 0E (3 ) 7 RS0=0. ...(14) 


The modified boundary conditions are 
f= 0:u=0 forally 
u=—= fr” for y= ‘| ...(15) 


r= 0's 
u= Ofory= 1 
and 
t= 0:6= 0 forally 
eee = | ‘ 
z = 0 fory = 0 f ...(16) 
P>20> 
6=Ofory=1. J 


3. MeTHOD oF SOLUTION 


It is difficult to solve equations (13) and (14) which are of order 3. However, 


one of the possibility is to derive one solution for small values of Rc (< 1) and then 
we assume 


u= Up + Rou, + O (R?). etl 


Substituting (17) in (13) and equating the like powers of Re we get 
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Ou eu a7 u 

tah fo oan (areal ate _ Se 

ot ze cy oy? ...(18) 
Cu, ( Ou, Oto Cu, Fu 

—— > tee, pao fs 0 a 

or - dy? + dys ) ay? eT ae «n{19) 


The boundary conditions for eqns. (18) and (19) are obtained from eqns. (10) and (17) 
as follows: 


t=0:u= 0, u,=0 ¥ y 
iu, = Ofory = | ...(20) 


Uo 
Yess ol © 
Uu = Ou, = Ofory = 1 


Equations (18), (19) and (14) subject to the conditions (20) and (11) are solved using 
Galarkin technique. For this the following infinite forms for uo, u, and @ are pro- 
posed which satisfy the initial and boundary conditions. 


Uy = t™(1—y) + aty(l — y) +4, ty? (1 — y*) +... 


az t® y3 (1 — y)§ +... as(21) 
uewbhyl—-ytheyv(l —y? +b, 0y2(l—y)e +... (22) 
6yect(l —y)+et y(l — y?)? +, ey? (1 — y*)® +... ...(23) 


where a, 5; and c; (i = 1, 2, 3) are arbitrary constants to be determined. Substituting 
(21) — (23) into eqns. (18), (19) and (14) and neglectings a, b; andc; fori > 4, the 
defect functions D, » D,, and De are obtained which are minimized by the technique 


of orthogonalisation and result in the following nine double integrals 


if pl 
ff Dy tiyi (1 — y)i'dtdy = 0 
00 


i Wows! 
J J D, tiy'(l — y)idtdy = 0, j=1,2,3 ...(24) 
0 1 


| 
| 
> 
| 
es : 
f {Det ji (1 — yt) dtdy =O. | 
60.6 («0 
J 
It is note worthy here that ¢ € [0, 1]. 


Performing the integrations, we arrive at the nine algebraic equations involving 
the parameters a), b, and c,(j = 1, 2,3). These nine linear equations are solved 
which give the constants a,, b) and c, and hence the velocity field « = uy + RK. u, and 
the temperature field 0. 


The non-dimensional shear stress +x» is given by 


SS 5 Ce | a fu , Ou } 25 
aa, AT* = = — Re eee aah vae(25) 





git eee 
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The shear stresses at the plates t) = txy | y.o and t1 = Txy | yo1 are calculated. The 
rate of heat transfer at the plates 
6 060 
Nug =~ SS and Nip te ae 


Yun 


are also calculated. 
Yum) 


4, RESULTS AND DISCUSSIONS 


The flow characteristics are studied for two positive values of 7, 1. €., 
(i) n = 1, constant acceleration 


(ii) m = 4, variable acceleration. 


It is observed from the Fig. | that the flow field exhibts almost a diametrically opposite 
behaviour for values of R. greater than zero. It is further observed that the velocity 


1. 0-05 5.0 0.0 
1, 0.05 50 0.9 
tm. 0.05 5.0 0.” 
V- 0.15 5.0 0.10 
Vv. 0.15 10.0 9.10 





Fic. 1. Effect of R and R, on velocity field, 
n= 1.0 


Fic. 2. Effect of t and n on velocity field, 
R= 5.0, Re = 0.1 


rises sharply with increasing values of f and this rise is quantitatively presented as 
152%. Torecord the relative influence of the two parameters ft and n on the flow 


pattern, it can be said from Fig. 2 that time plays a more effective role in increasing 
the velocity than n. 


It is noticed from Fig. 3 that an increase in R, is to decrease the temperature at 
all points and also an increase in time decreases it further. 


Figure 4 shows that an increase in R is to increase the temperature at all points. 
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t Re 

I. 0.05 0.0 

tl- 6.05 0.05 


il. 0.05 0.10 
IV. 0.15 0-10 a 







4.0 








3.0 
e 
2.0 
1.0 
9.0 2 
0.2 C4 0.5 2.5 1.0 
4 ay 
Fic. 3. Effect of R, on temperature field, Fic. 4. Effect of R on temperature field, 
R = 5.0,n = 1.0, S = 0.1, co = 0.0, R, = 0.05,” = 1.0, S = 0.1 
E= 0.01 E = 0.001,o0 = 0.5 


Figure 5 shows that the temperature field assumes negative values for S = 1.0 
and S = 0.5. Further, it is observed that the temperature decreases with the decrease 
of source strength from 1.0 to 0.5 but for S = 0.1 it increases sharply. 


Figure 6 shows thatthe shearing stress at the lower wall decreases with an 
increase in time for both variable and constant acceleration. It is also noticed that 
t,, in case of constant acceleration is greater than the variable acceleration for fixed 





etiam RE 
* =—0.4 


Fic. 5. Effect of S on temperature field, R = 10.0, Fic, 6. Shearing stress at the lower 
t = 0.05, R, = 0.1, n = 1.0,0 = 0,2, E= 0.01 wall, R = 0.1 
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value of t. Whereas the shearing stress at the upper wall +, depicts an opposite 
behaviour for both ft and n (Fig. 7). Further it is observed that the elasticity of the 


liquid increases both t, and *,. 


0,05 n=05 


t= 

It. 33,0.) n=0.5 
Il t= 0,05 n= 1.0 
TV) -t =.'0) n= 1.0 





Fic. 7. Shearirg stress at the upper wall, R = 0.5. 


TABLE [| 
Values of the rates of heat transfer, Nug and Nu, S = 0.1,0 = 0.1, E = 0.01 








n 0.5 1.0 
R t R,/n Nuo Nu, Nuo Nu, 
0.05 0.0 0.00004 0.00160 0.00003 0.00131 
0.05 0.00006 0.00194 0.00006 0.00193 
0.10 0.00004 0.00159 0.00003 0.00129 
0.5 
0.0 0.00019 0.00321 0.00375 0.00262 
0.10 0.05 0.00024 0.00389 0.00024 0.00387 
0.10 0.00019 0.00319 0.00015 0.00258 
0.0 0.00015 0.00628 0.00012 0.00520 
5.0 0.05 0.05 —0,00112 —0.01647 —0.00095 —0.01449 
0.10 —0,00047 —0.01097 —0.00032 —0.00737 
TABLE If 


Showing the effect of S on the rates of heat transfer, Nuy and Nu, t = 0.05, Re = Oe. 
n= 1.0,c = 0.2E = 0.01 


EE... 





R S Nupo Nu, 
0.1 0.00004 0.00140 
0.5 0.5 0.00004 0.00140 
1.0 0.00004 0.00141 
0.1 0.00014 0.43116 
10.5 0.5 0.00735 —0.65259 
1.0 0.00496 —0.13865 


eee 
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Table I describes the effects of R, Re, tand on the Nusselt number. Nu, and 
Nu, both increase when R; takes the values from 0.0 to 0.05. The rate of heat transfer 
falls for Re = 0.1. From Table II it is ascertained that the increase in source strength 
maintains almost a constant rate of heat transfer at both the walls for R = 0.5. 
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ON RAYLEIGH WAVES IN GREEN-LINDSAY’S MODEL OF GENERALIZED 
THERMOELASTIC MEDIA 


N. C. Dawn AND S. K. CHAKRABORTY 
Department of Mathematics, University of Burdwan, Burdwan 713104 


(Received 24 August 1987; after revision 23 March 1988) 


Rayleigh wave formulation in a half-space exhibiting generalized thermo- 
elasticity properties of Green-Lindsay model has been made. The case of 
Lord and Shulman’s generalized thermoelasticity is shown to occur asa 
special case. Approximations to the frequency equation have been made for 
different ranges of the parameters involved. In the case of small coupling, 
a numerical study has been made to show the effect of the relaxation para- 
meters on the amplitude attenuation factor and also on the velocity of 
propagation 


INTRODUCTION 


Thermo-elastic Rayleigh waves in semi-infinite isotropic solids have been studied 
by Deresiewicz', Lockett*, Chadwick and Sneddon*, Chadwick‘. Chakraborty and 
Pal® extended the problem to transversely isotropic medium. These surface waves are 
found to propagate with a speed in general dependent on the thermal properties as 
well as on the wave-length of the waves. In the case of small reduced frequency, how- 
ever, the thermal effect on the wave speed is found to be negligible. 


The physical foundations of the above problems were the formulations of Biot® 
and Lessen’. The coupled heat conduction equations used by them permitted infinite 
speeds of propagation of thermal signals. Lord and Shulman’, proposed generalized 
thermo-elasticity equations, in which they modified Fourier’s heat conduction equa- 
tion by taking into account the time needed for acceleration of the heat flow. This 
leads to changes in the energy equation of a coupled theory of theoremo-elasticity by 
the incorporation of a relaxation term. Green and Lindsay® proposed generalized con- 
stitutive equations starting from the thermo-dynamical principles without altering 
Fourier’s heat conduction equation. The Green-Lindsay thermoelasticity equations 
involve two new parameters. Both the Lord-Shulman and Green-Lindsay formulations 
give finite speeds of thermal wave propagation and also a “second sound effect”. How- 
ever, Green-Lindsay formulation appears to be theoretically more satisfying and Lord- 
Shulman case follows asa special case of Green-Lindsay in many problems. In the 
linearized case, according to the Green-Lindsay theory, the relaxation parameters a and 
a* in the stress-strain temperature relation and the equation of heat conduction 
respectively are independent but the inequality a > a* > 0 is to be satisfied for 


RAYLBIGH WAVES 277 


uniqueness of temperature and the speed of heat conduction which depends on a* is 
finite only if the stresses depend on the time rate of temperature. 


Nayfeh and Nemat-Nasser"! used the Lord-Shulman theory to study plane thermo- 
elastic surface waves in a half space. On the basis of Green-Lindsay model, Agarwal’? 
studied plane thermo-elastic waves—their propagation and stability. The two results in 
general are different. They coincide however, in the special case when the two relaxa- 
tion parameters are equal. 


The present paper is concerned with the problem of thermoelastic Rayleigh 
wave in semi-infinite solids of Green-Lindsay’s model. It is found that the relaxation 
parameters a, a* for Green-Lindsay theory contribute to the frequency equation terms 
of order higher than X-"/*, % being the reduced frequency. It is seen that when the two 
relaxation parameters are equal, the results reduce to the case studied by Nayfeh and 
Nemat-Nasser™ based on Lord-Shulman’s theory. 


Approximations to the frequency equations for different ranges of the parameters 
involved have been obtained. A numerical study of the change of the amplitude atte- 
nuation factor and also of velocity of propagation with the relaxation parameter a have 
been made in the case of small thermo-elastic coupling. 


Basic EQUATIONS 


The equations governing linear thermoelastic interactions in a homogeneous and 
isotropic solid free from body force and heat sources, as proposed by Green and 
Lindsay® and Green” are : 


(a) the strain-displacement relations : 
2e;, = ui,j + Uy eur) 


(b) the stress-strain temperature relations : 


ty =AAS,+2pey — y¥(T + aT) 5, vetle) 


(c) the equation of motion : 


+> — >> ° > 
pVutAt+pyvdivu-—yv(T+aT)=Ppu ia) 


(d) the equation of heat conduction : 


«VT — pe, (T+ a* T) —yTMdiv u=0 (4) 


where 
a ao = 0 


e;, = cartesian components of the linear strain-tensor 


ty = cartesian components of the linear stress-tensor 
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5;, = Kronecker’s delta 
A = u,,; = dilatation, 
y = (3A + 2) a 3A, » = Lame’s constant 
a, = coefficient of linear thermal expansion 


« = coefficient of thermal conductivity, 


P = constant mass density 
cy = specific heat at constant volume 
T = 


the change in the absolute basic temperature 7). 


ate 


| 


thermal relaxation times (constitutive coefficients). 


FORMULATION OF THE PROBLEM 


For Rayleigh type waves in the half space z > 0, using the representation of 


displacement components : 

a eae ae 

me eee, 

where ¢ and y are functions of x, z and f, eqns. (3) and (4) are satisfied if 


ed ad o2 v . 
Fm ch reel ye Oe) 


ey ¢ ary au 
af. oS eR r 








and 











A 2 2 P 
Pe(T+at +7 5 -| et sar bae( Sa 


with c) = (A + 2y)/P, cl = p/P. : 


Using the quantities 





ae = x](¢,/@*), 48 = z)(c,/*), t' es t w*, 
‘ Ux ‘ uz , ae T 
Fe => ———. = 
om (Cjaty' safety ; 
pc? 
, d , wv , 
= Gare case 14 Saar, ah se ia® at 


where w* = 


equations in dimensionless form : 


<2 {5) 


...(6) 


te) 


oT 
+ d z* ) 


.. (8) 


Pc, Cc /<, in eqns. (6), (7), (8) and suppressing the primes, were obtain the 
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MS Ox ey 2) ag , ad 
ne oN pet gs, tay pot 
a a? db ad 
Or = (F5+ St) - tan ...(10) 
; = a a o a: oT a2 T 
ih > ——4 | —— — = ene 
a ee a call ax? + ey ) axt + 9a St 
and 
o2 wv 1 e2 v eo y 
ato (45 — ) (12) 
where the thermo-elastic coupling is given by 
2 a 
see To and vy? = 
HP IA or : 
gis e “s ate 
Stress and thermal conditions of the problem on the boundary z = O are tix = tz 
Ped sndeee fo 0: (14) 
0z 
SOLUTION 


For thermoelastic surface waves in the half-space propagating in x-direction the 
functions {7, ¢, ¥} may be taken in the form 


{T,¢,v¥} = {T (2), b (z), b (2)} exp i(n x — X 2). act 15) 


A AN A 
Substituting (15) in (9) — (11) and remembering that 7, ¢, ) ~ 0 as z— c¢ for surface 
waves, the solution is obtained as 


¢ = [A exp (— 7 B, z) + Bexp(—7 B: z)]. expi(ax — X ft) .»(16) 

~ = C.exp([—n8;z+i(4x —%2)] 2612) 
and 

T= er [A {X? + »? (8; — 1)} exp (— 1 8: 2) 


+ B{x? + 4° (85 — 1)} exp (— 1 B, z)]. expi(gx — X14) ...(18) 


where B; = (1 —c* v’),c? = x?/n? and B,, B, are those roots of the following equation 


for which Re (B) > 0: 
a6+ B[—2+e(itat+ea+i(l +OMI+1 —c? (A+e" 


+ea)+ _ (l+6}t+et@* + —)=0. ...(19) 
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On using the boundary conditions (14), A: B: C are eliminated giving rise to the fre- 
quency equation 


(2 — v?c*)? (8B + PB, +81 B. — 1 + c*) — 4B, B: Bs (8, + 8») 


SP o ma + B,) (2 — v2 c?)? — 4 B3 (8, B, + 1 — c*)}. 
...(20) 


DisCUSSION 


Equation (20) along with (19) determines the velocity of thermoelastic Rayleigh 
waves in a relaxing medium. In order to have an idea of the effect of the relaxation 
parameters on the velocity of propagation we shall study below in detail the special 
case of small thermoelastic coupling : 


Case (a): Small Thermal Coupling 


For most of the materials ¢ is small at normal temperature. Hence we may make 
an approximation of the frequency equation assuming « < |. 


For « < 1, we get from (19) 





i 
oa 
B= (1 eye [ 1 — 5 ( ik ) 
Qa -—>)0-e) (21) 
swt bee bog 
(1 — a*—i/X) {1—c? (a*+/]% / 


The frequency equation (20) with h = 0 thus reduces to an equation in which 
the parameters a, a* are involved : 


aH —eat + Syn 4 t 


Q—vey[l—c(at teat FU + 9 +0 -eF le (a 


| yup gq — _f c* (a + i/X) 2 \2)1 
FR ae) {ioe ee ee 


]—c¢?)}/25] —¢2 E vip wha c? (a + i]X) 
(1—c*)'/*{1 ~—c¢ (a* + 9 ie {1 2 ieee HO) 


—p2)l /2 / vit 5 X 
[(1—e2)'2 {1 —e ee + (1-8 (a*tifppe 


Toe. c? a* + i/X) (a + i/X) x 
{1 + /2 Wary lap! = O- ...(22) 
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If we put 


C= c¥? + ¢ (Sy + 76,) 
where c* is the classical Rayleigh wave velocity and &; and =, are two reals depending 


on the reduced frequency % and a, a*, 


then 
a x ene eT 


c* 2c#20— Dc? 











The velocity of propagation being (c* + ais ) the waves exhibit dispersion and 


the amplitude-attenuation factor = exp [ add a ] with —& < 0. 








2c*8 


Substituting the above value of c? in (22) and solving numerically for different 
values of the parameters we have calculated the velocity of propagation and the ampli- 


tude-attenuation factor. In Fig. 1, we have plotted velocity of propagation (X) 10? 
against the relaxation parameter a for « = 0.05, (4) = 0.1 and a* = 0.2 and in Fig. 2 


amplitude attenuation factor (X) 10* has been plotted against the relaxation parameter 
a for the same value of e, % and a* for the case x = 1. 


From Fig. 1, it is seen that the velocity of propagation decreases gradually as the 
relaxation parameter a increases, though the change is very small. This may be inter- 
preted as solving down of Rayleigh waves with increasing heat conductivity (Nayfeh 
and Nemat-Nasser"’, p. 55). 


Similarly, in the Fig. 2, it is seen that amplitude-attenuation factor decreases with 
increase in the relaxation parameter a. 
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Amplitude attenuation factor x10° 


2 03 .04 65 0607 -08 TO 


Qa —-_—_ 


Fic. 2. 


Case (b): Small reduced frequency % < }. 


It is seen that the characteristic frequency «* is much greater than the frequencies 
attainable even in experiments employing ultrasonic pulses. Thus we may assume 
o {o*,i.e. %< 1. Expanding the terms of the frequency equation in powers of 
the reduced frequency %, the frequency equation (20) can be written as 





. 2 an : 
Qa rioy ic = (i+ ) toa 140) ok eaee 
; (1 +i) (1 —c?(1 + at + ea) +. a* c4} 
(a* + €a) + 7a/d /2 TT Pe Cee ee ee ud 


— 4(1l- eyr| — fe (1 +e)P Ud +e —c*)P2 
(451) cs PS ele 
+t V/2 x12" (1 + et? Heat Lise) 





v, {3 — 2c? (1 + a*+ ea) + a* ct} — c? {2 — c? (1 + a* + ca)} 
, 2(1 + ei (1 ao a ‘ed hte 


4 Q+i) “4 {1 —c*(1 + a* + ea) + a*c4} 
2/2 Sake (GI ge 9 bbs ] 





h 1—j c (1 + « — cP (1 + i) 
=. — = ——_ , — /2 ee eee 
i Jo) eee G+eor t+ 2y2 


Ye yee a * 
ee eee } (2 — v2c2)? — 4 (1 — y2c212 


(equation continued on p. 283) 


RAYLEIGH WAVES 283 





I — | £ x ] 1] xy /e 
x | OL yi . (1 a €— oye a D oe ae 
Cp a i LIME we es oA) Sh : 
{aoe + 1 e*}}, .-6(24) 


From (24), we therefore find that the relaxation parameter’s ‘a’, a* are involved 


only in the terms of order higher than X-2._ Hence the effect of these parameters is 
negligible for small X. 


The results corresponding to zero flux of heat across the boundary may be obta- 
ined by putting A = 0. Therefore putting A = 0 in the frequency equation (24) and 
retaining terms only O (%~") we get 





- l1—e-—c¢ 
Piss 28 WEE ey Peete 12 \ 1/2 
( ¥¥ ¢3) 4 (1 — v? c?)'/?, Terai ae Be eR 
Equation (25), is identital with the corresponding equation in coupled thermo- 
elastic medium’. 


Case (c): Lord-Shulman Case 


It is to be noted that eliminating the temperature from the displacement equation 
(3) and the heat conduction equation (4) and then putting a = a* results in an equa- 
tion for displacement which is the same as in Lord-Shulman’s case with ‘a’ as the 
relaxation parameter. Therefore in Green-Lindsay’s theory, all wave problems and 
source problems reduce to Lord-Shulman’s case with a = a* playing the role of re- 
laxation parameter. 


If we put a = a* then the characteristic equation (19) and the frequency equation 
(20) coincide with the case corresponding to Lord-Shulman theory’?. 
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Kopal!’* obtained expliciting expressions for Roche harmonics associated with 
the stellar models of a star distorted by rotation (solid body) forces Follow- 
ing Kopal’s approach we obtain explicit expressions for Roche harmonics 
associated with the stellar model of a star distorted by the differential rota- 
tion (the law of differential rotation has been assumed of the form 
© = 5; + b, s°, where is the angular velocity of rotation of a fluid element 
distant s from the axis of rotation and b,, 6, are certain constants). 


1. INTRODUCTION 


Kopal’ introduced a family of new auxiliary functions, which he called the 
Roche harmonics, generated by the solution of the Laplace equation 
ve = 0 kok G5 
derived in terms of the Roche curvilinear coordinates, which are associated with the 
Roche equipotential surfaces in the same way as the spherical harmonics are associ- 
ated with a sphere. Roach® discussed some more theoretical aspects of Roche 
harmonics. In fact, as was pointed out by Kopal (op. cit.), spherical harmonics may 
be regarded as a limiting case of the more general Roche harmonics, and the latter 
may be more appropriate for the study of many problems arising in double star 
astronomy. Kopal obtained explicit expressions for Roche harmonics for rotatio- 
nally or tidally distorted stellar models. In the case~of Roche harmonics associated 
with equipotential surfaces which are being distorted by the rotational forces, Kopal 
obtained explicit expressions by taking into account the effects of solid body rotation 
only. Singh*, Mohan, Singh® extended the analysis of Roche co 
into account the effect of differential rotation. 
been assumed to be of the form w» = b, + bs? 
an element distant s from the axis of rotation, b, 


ordinates taking 

The law of differential rotation has 
» Where w is the angular velocity of 
and 5b, are constants. 


In this paper we present the analysis of Roche harmonics for stellar models 
distorted by differentia] rotation, using the above cited law of differential rotation as 
« = b, + b, s?, to obtain the explicit expressions of Roche harmonics associated with 


Roche equipotential surfaces distorted by differential rotation. 


This is performed in 
Section 4. 
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For completeness the work of Kopal! is respected in Section 2 considering 
solid body rotation effects only. The system of Roche coordinates for the Roche-model 
distorted by differential rotation as obtaine dby Mohan and Singh® are presented in 
section 3. Concluding remarks are given in section 5. 


2. Expxicir ExpreSssIONS OF ROCHE HARMONICS FOR ROCHE EQUIPOTENTIAL 
SURFACES DISTORTED SOLID BY Bopy ROTATION ForRCES 


Following Kopal' it can be shown that, 
—g= : + nr? (1 — v*) = constant »+(2) 


represents the equipotential surface of a star distorted by rotation (solid body) forces. 
Here n = @?/2, » being the angular velocity of rotation (solid body rotation), r is the 
dimensionless measurement of distance and v = Z/r is the direction cosine with respect 
to Z-axis, the axis of rotation. The surfaces generated by setting & = constant in 
eqn. (2) are referred to as Roche equipotentials. 


Choosing & as defined in eqn. (2) as one coordinate Kopal has shown that the 
second Roche coordinate % is given by 


n= ¢. (3) 


The expression for the third Roche coordinate ¢ is not possible in closed analytic 
form. Kopal obtained expression for this third Roche coordinate ¢ in ascending 
powers of the angular velocity of rotation as 


cos f = z (2n)’. r3) X; (v) ...(4) 
j=0 
where X, (v) = 1, X, (vy) = — 41 — v?), while for j > 1 all subsequent Y, (v)’s can 


be generated with the help of the recursion formula 
37 X, + (1 — v4) (Xa) — 3G — 1) Xj] = 0. 
Here prime denotes differentiation with respect to v. 


Kopal also obtained expressions for the metric coefficients 4,, 4, and hs in terms 
of polar spherical coordinates upto second order terms in square of angular velocity 
of solid body rotation. 


These are : 
h, (é, 0) = ré (1 + 4nré sin’s — 1 n? rf sin’ (22 — 85 sin’r) +-.-] 
h, (é, 0) = ro sing [1 — 4 nr (2 — Ssin’s) + $n? ro (2 — 50 sin’ 
+ 7 sin‘s) +...; 
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hs (€, 6) = ro[1—} nr8 (2 — 7 sin’%) + $n? r§ (2 — 88 sin? C 
+ 145 sin't) + ...]. ...(5) 


Here r, = 1/— denotes the mean fractional radius of the equipotential surface 
= constant. 


Making use of /,, h, and /s as given by (5), Kopal? has shown that the Laplace 


equation (1) becomes 


; 2 
vie= [l—an (l- %) +3 my (=) 6 Fe loy I 

6 2 Op — Od Z 

or, (re $E-) — anes + [1 — > me (5 — 7v, ) 


l o2 
— gy (43 — 194 vy + 143 | 4) ? 





0 


9 
_|, 4 3. cee 2 4 
+ 3 My Ys tener, ioe tie ees 





9 
6g 
2vo eB on 0 (6) 
where 
v = cos &, ro = 1/E. 
Assuming a series of the form 
¢=Za,r/ R 
E a, wen, on od 
where 
R, = P,(v)+ nry XY (vo) +n? cr XY (,) +... -3(8) 


Kopal has shown that for j > 1, gs (vo) and > a (vo) are found to be of the form 


X” (v0) = — (1 — v3) (1 — 5v8) ..(9) 
X~ (vp) = — 4(1 — v3) (11 — 25 v8) vp .. (10) 
l 
xX) = % (1 = v6) (9 — 120 vg + 175 v8) sat hi) 
x (y,) => (I 2 2 
2 (%) == (1 — vo) (17 — 98¥8 + 105v4) vo ...(12) 


( ] 
Xo (yg) 6 (1 — vo) (21 — 172v8 + 175 v4) ...(13) 
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etc. 


Equation (8) with (9)-(13) constitute the explicit form of the Roche harmonics associ- 
ated with Roche equipotential surfaces (2) distorted by solid body rotation forces. 


3. ROCHE COORDINATES FOR ROCHE MopeL DisToRTED BY 
DIFFERENTIAL ROTATION 


For the Roche model of mass M, rotating according to the law 
w=b + d,s° .. (14) 
the equation of hydrostatic equilibrium may be written in the form 
dQ = dy + 3 w* d (s*) me AI) 


where Q denotes potential at a point P distant r from the centre of the 
star. v = GMJr, is the gravitational potential, » is the angular velocity of rotation of 
an element of the fluid distant s from the axis of rotation. 


On substituting (14) in (15) we have 
dQ = dv + } (b? + 2b, b, s* + b3 54) d (s*). 
On integration and simplification this gives 


= 2M sg (xt + yt) [OE + by by (x + y*) +4 :0E OP + y°F. 18) 


F 


In spherical polar coordinates 


x =rcos¢sin& = ra 
y=rsin¢gsin@ = rp att) 
z =rcos 6 =v 
Equation (16) may be expressed in non-dimensional form as 
p= Leg —wbtt b.o,(1—v) + bor — 94 - 


where & -= RQ/GM is a non-dimensional parameter denoting potential and is NOW 
non-dimensional angular velocity in units of GM/R?, R being the unit of distance. 


The surfaces generated by setting & = constant in (18) are referred as the Roche 
equipotentials. 


Now if we take ro = 1/& as our first approximation to the distance of the equi- 
potential surface from the centre, Mohan and Singh’ have shown that 


rero(1 t+ 3r3 (1 — v5) (65 + A: 0 + w$)] (19) 
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where 
vo = cos Cand w = db; + bz rf (1 — vo). 


In the system of Roche coordinates (&,7,¢), the €-coordinate is defined by Roche 
equipotential surfaces of the form (18) while coordinates 4 and ¢ are defined by the 
requirement that they are orthogonal to € as well as to each other. In this triple 
orthogonal system of Roche coordinates, Mohan and Singh® have shown that the 
second and third coordinates are given by 


A 


ne yt ..-(20) 


and 


Tae 
costa y[1— CTI" (sue + 120+ 86) +. ] 21) 


whereas the expression for n given in (20) is exact, the expression for ¢ obtained in 
(21) is correct upto second order terms in angular velocity ©. 


Mohan and Singh’® also obtained the explicit expressions for the metric coeffici- 


ents /;, h, and hs correct upto sectond order terms in », which are found to be 


hs (Ete [1+ one 8 4 by tip + BE) + .. 4 


h, (E,5) = rq sinG [1 - (15 w§ + 126, mw) + 8b?) 


70 


r’ sin’ 


+ ThSIDY (65u8 + 9b: %» + 51b%) + =| 








3 
hs (&, ©) = ro| 1 — 2% (150 3 — 48 b, w» — 3264) 





— sin®t (215w3 + 11, @» + 1963) + | (22) 


4. ExpLicir Expressions oF RocHe HARMONICS FOR ROCHE EQUIPOTENTIAL 
SURFACES DisTORTED BY DIFFERENTIAL ROTATION FORCES 


Using expressions of /y, A, and hs as given by (22), the Laplacian 


.2 ] | ae (i Fa) h,h, @ afi haa 
Vo he Laie) + ay a) tae) 


now becomes 








9 ro 


y= [1 —_ 105 (1 — vo) (2805 + 70 by wo +.70 1) | a(n g ) 


r 
° aro 


(equation continued on p. 289) 
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arg $ 
. ie (9003 — 12b, &) — 852) 


are 2ro bs 


+ SR? (1 = v3) (18005 — 1264) | 


ro ss 
- [2 + a (60005 — 192 b\w) — 128 b}) v2 
+ (1 — vG) b, r§ ((520 (1 — vg) 
+ 480) w. + (236 (1 — v3) — 144) bj) ] Yo ie 
‘ I = mas — v2 (43002 + 22 b, wy + 3862) 
+ 130 w2 + 1185, we + 102 3 i _ eae 
Therefore the Laplace equation (1) takes the explicit form in this case as 
ji- a (1 — v8) (280 2 + 705, wy + 706%) | 5— ( 73 ~ ) 
105 ero 


2re:0 


2 
re (9005 — 1261 %» — 8bi) + p5~ (I — va) 


105 





ad 
(180 = 12b,) | rs Or, 
8 

as | 2 + FB (600 of — 192 b, w — 12864) vi 

+ (1 — v8) by r2 (520 (1 — v3) 

4. 480) @, +236 (1 — v8)-— 144) b,) | ; 

4 E - a — v3 (43002 + 22b, w + 3863) 
Ee (23) 

4 130 3 +1185, w, +1023 ] (v8) GF = 0. 7 


As in section 2, the solution of this equation can be obtained in series form where 


¢=Zajr R; 
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and setting 

R, = P;(%) + aD reek egy cee (24) 
It can be shown (consistent with the adopted scheme of approximation) that the 
functions x (vo) assume the forms because of the influence of the terms pertaining 


b3 and b, bs, as given by 


X® (v9) = — 1.8571 vé + 2.1905 v2 — .3333 ae) 
X$” (vo) = — 4.6428 v8 + 6.5714 v8 — 1.9286 vo ...(26) 
X! (v9) = — 10.8333 v8 + 17,9762 vé — 7.6428 v8 + .5 20 
X{) (v9) = — 24.375 v2 + 46.458 v8 — 25.7440 v3 + 3,6607 vo; ...(28) 
Dhl C9 ey im a Be et ...(29) 
Xs” (v9) = — 6.75 vB + 9.6 v3 — 2,85 vo .. (30) 
X5" (vo) = — 15.75 v8 + 26.25 vé — 11.25 v2 + .75 me) 
X,” (vo) = — .33.4375 v3 + 67.8125 v8 — 37.8125 v8 + 5.4375 vo...(32) 


etc. 


Equation (24) with (25) - (32) constitutes the explicit form of the Roche harmo- 
nics associated with the Roche equipotentials surfaces (18) up to second order appro- 
ximation of the differential rotation terms, which indeed, is the outcome of our 


analysis. 
5. CONCLUDING REMARKS 


The influence of b? terms on xe (vo) is same as obtained by Kopal*. In the 
present analysis we could only present the results which are because of the influence 
of the terms pertaining b3 and b, b,. This is because of Our approximation scheme 
for angular velocity » upto second order. The consideration of higher order terms 


in ©, should lead to the appropriate formulation of this problem which we intend to 
investigate this problem in subsequent study. 


It may be pointed out that although we have studied here the problem of Roche 
harmonics associated with equipotential surfaces by assuming the Roche model of the 
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star, the present method of Roche coordinates can be used also when some more 
realistic structure is assumed for the interior of the model. We can still either 
approximate the distorted equipotentials of such stars by Roche model or use their 
more realistic forms by using the system of Clairaut’s coordinates (cf. Kopal®’’), when 
the law of differential rotation may be assumed of the form o = b, + b, s°. 
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The two-dimensional problem of incoming waves against a vertical cliff 
ina two-fluid medium is considered in this paper. The velocity potentials 
describing the motion in each of the fluids are obtained. Known results 


for the one-fluid medium are recovered when the density of the upper fluid 
is made to vanish. 


1. INTRODUCTION 


Stoker*”* considered the problem of incoming surface water waves in a deep ocean 
bounded on one side by a vertical cliff. He solved the two-dimensional problem by 
employing a technique based on the complex variable theory and the three-dimensional 
problem by a reduction prccedure. No reflection of waves by the cliff is assumed so 
that a source/sink type tehaviour of the potential function there is necessary to account 
for this and as such the wave amplitude is required to be logarithmically infinite at the 
shore line. However, if the effect of surface tension at the free-surface is taken into 
account, this is not necesary, so that the potential remains regular there. Packham! 


solved the corresponding two-dimensional problem by using a technique based on the 
Fourier sine transform. 


The present paper is concerned with the two-dimensional problem of incoming 
internal waves at the interface of a two-fluid medium against a vertical cliff in the pre- 
sence of interfacial tension thereby generalising Packham’s! problem. The velocity 
potentials in both the fluids are obtained in a straight forward manner. 


2. STATEMENT AND SOLUTION OF THE PROBLEM 


We consider the irrotaticral motion of two inviscid incompressible fluids of in- 
finite horizontal extent towards the right and bounded by a vertical cliff on the left and 
of densities P, and p. (< P,) respectively under the action of gravity only. The lower 
fluid extends infinitely downwards while the upper extends infinitely upwards. We 
choose a rectangular cartesian co-ordinate system such that the y-axis points vertically 
downwards into the lower fluid, then the cliff is the plane x = Oandy=0,x>0 is 


WATER WAVES 293 


the position of the undisturbed interface. The velocity potentials ®, (x, y, t) and 
”, (x, y, t) for the lower and upper fluids respectively, under assumption of time-har- 
monic motion, can be described by 


©; (x, ¥, t) = Re [9, (x, y exp) (— ist)](j = 1, 2) 
where ¢ is the circular frequency. Then 9; s satisfy 


V* 9, = O in the respective region fe) 


Piy Pry 
Kei + 9 ~ 5 (Ker +9n) + MI pita = 0Oony=0, x>0 


Poyry 
ct22) 
gx = Oonx = 0(y > 0 forj = 1 andy < 0 for; = 2) ey 
p, S remain finite at the origin .. (2.4) 
V9, > Oasy > + 0 mil) 
(the upper sign is for j = | and the lower sign for j = 2) and finally 
9, ~ + exp(F ky y — ikyx) as x > 00 snt2,0) 


where ky is the unique positive real root of 
k(li+ M'k?)—-L=0 test) 
M’=M/(1—s) L=K(1 + s)[(l —s), K = o]g 


s = P,/P,, M = T/(P, g). 
T being the interfacial tension and g the gravity. 

Here, (2.1) is the equation of continuity in either of the fluids, (2.2) is the line- 
arized kinematicel and dynamical conditions at the interface, (2.3) is the condition at 
the cliff, (2.4) is the condition at the shore-line, (2.5) is the condition of no motion at 
infinite depth and height, and (2.6) is due to the incoming nature of the waves as 
x — co moving towards the cliff. 

To solve for 9; (x, ¥) we write 
9, (x, y) = + 2 exp (F ko y) cos (Kyx) + wy; (x, y). Race, 
Then ¥,’s satisfy the same eqns (2.1) to (2.5) as satisfied by p,'s and 


b, (x, Y) > F exp (F ky y + ik, x) as x > 0. em RS) 


(2.9) states that /,’s behave as outgoing waves as x > co. 
Solution of ,’s satisfying (2.1) to (2.5) is given by 


%, (xy) = +e H eee - cos (kx) dk (j = 1,2) (2.10) 
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where the path of integration is indented below the pole at k = ky to account for the 
outgoing nature of /;’s as x > co and c is a constant to be chosen such that (2.9) is 
satisfied. We may note that ¢,’s given by (2.10) remain finite as (x* + y*)/* +0 so 
long as T > 0. However, for T = 0, (2.10) exhibit a logarithmic singularity as r > 0 
(cf. Yu and Ursell*) which accounts for a source/sink type behaviour at the shore-line 
in the absence of interfacial tension as stated in the introduction. 


Alternative representation for ); (x, }’)’s are given by 


ni : 
v,(x y= + ¢| EE exp (= ky y + ik x) 


00 
+ | Sore ne (l—M k?) cos ky a ZL sin ky} dk I. 
0 


wale 
To satisfy (2.9), we must choose 
i eis 
c= a, + 3M’ k3). wasteste) 
Thus we obtain finally, 
D; (x, y, t) = + exp (F ko y) cos (ky x + ot) 
P co 
ap eee 3M’ k3) | 
0 
exp (—kx) {k (1—M’' k*) cos ky + L sin ky} wa(2.13) 


kK? (1 — M' k?)? + L? 
In the absence of interfacial tension (2.13) reduces to 
®, (x, y, t) = + exp(F Ly) cos (Lx + 6 £) 


, Sin ot [spe kx) 


ke {k cos ky = L sin ky} dk....(2.14) 


7 


Putting s = 0 in the expression for ®, in (2.14), the potential function obtained by 
Stoker?” is recovered. Again, putting s = 0 in ®, into (2.13) we obtain 


®; (x, y, 1) = exp (—ko y) cos (ky x + st) 


sin of f° 
+ 1+ 3Mi exes Tae o 
an KG) k2 (1 — Mk y+ K 


{k (1 — Mk’) cos ky — K sin iy} dk 
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which was derived by Packham! using a different approach. 


4, DISCUSSION 


Potential functions representing incoming progressive waves against a vertical 
cliff in a two-fluid medium are obtained. Putting s = 0, the result for a deep ocean 
is recovered, more specifically, ®, reduces to Stoker’s?’® result when the surface ten- 
sion is ignored and to Packham’s! result when this is included. Also the potential func- 
tions in a two-fluid medium in the absence of interfacial tension are derived by simply 
putting T= 0. This problem can be further generalised to include the cases where 
the lower fluid is of uniform finite depth and/or the upper fluid is bounded by a hori- 
zontal rigid lid or a free surface. 
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